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ABSTRACT 

The  scope  of  this  thesis  is  to  evaluate  the  use  of  the  MINI  RANGER. 
LORAN  (TD)  and  LORAN  (DISPLAY)  navigation  systems  in  order  to  support 
the  collection  of  current  profiles  by  an  Acoustic  Doppler  Current  Profiler  Re- 
corder (ADCP). 

A  theoretical  error  analysis  of  these  systems  is  undertaken  in  order  to  estab- 
lish the  minimum  error  limits  which  might  be  expected  when  averaging  current 
profiles  over  time  frames  of  up  to  30  minutes.  ,>.,   , 

Experimental  data  was  collected  with  all  of  these  systems  in  the  Monterey 
Bay  and  was  analysed,  the  results  being  presented  in  this  thesis.  In  addition  GPS 
data  was  also  collected  but  time  has  prohibited  its  analysis  and  subsequent 
inclusion. 

The  results  show  that  because  of  ship  fluctuations  in  course  and  speed  there 
is  no  statistically  significant  difference  between  the  navigation  systems  when  we 
average  the  data  over  times  of  about  25  -  30  minutes.  However,  they  also  show 
that  both  the  MINI  RANGER  and  Bottom  Tracking  with  the  ADCP  can 
produce  reasonable  results  in  as  little  as  three  minutes,  although  the  ADCP  re- 
sults are  clearlv  biased. 


in 


TABLE  OF  CONTENTS 

I.  INTRODUCTION   1 

II.  DATA  COLLECTION    2 

■i 

A.  CRUISE   2 

B.  TIME  SYNCHRONIZATION     2 

III.  ACOUSTIC  DOPPLER  CURRENT  PROFILERS     5 

A.  INTRODUCTION    5 

B.  PRINCIPLES  OF  OPERATION  OF  THE  ACOUSTIC  DOPPLER 
CURRENT  PROFILER     6 

C.  THREE  DIMENSIONAL  CURRENT  VELOCITY  VECTORS    ...8 

D.  DEPLOYMENT  AREA    9 

E.  BOTTOM  TRACKING  OPERATION    10 

F.  VELOCITY  PROFILES    11 

G.  ACOUSTIC  DOPPLER  CURRENT  PROFILER  DATA  PROC- 
ESSING        11 

IV.  THEORY  OF  LORAN   C    13 

A.  INTRODUCTION    13 

B.  SYSTEM  DESCRIPTION     13 

C.  PULSES-PHASE  AND  CODES-CYCLE  SELECTION     13 

D.  TIME  DIFFERENCE  EQUATIONS    16 

E.  HYPERBOLIC  GRADIENT,  CROSSING  ANGLES    18 

F.  ERROR  ANALYSIS  FOR  LORAN  C     20 

V.  THEORY  OF  MINI  RANGER     30 

A.  INTRODUCTION    30 

B.  REFERENCE  STATIONS    30 


IV 


C.  RANGE  POSITIONING     30 

D.  ERROR  ANALYSIS  FOR  MINI  RANGER     31 

VI.  RESULTS    36 

VII.  CONCLUSIONS  AND  RECOMMENDATIONS     46 

APPENDIX  A.     STATISTICAL  ANALYSIS  FOR  THE  LC  408  DATA  47 

APPENDIX  B.     TABLES  FOR  STANDARD  DEVIATIONS  AND 
COVARIANCES  OF  VELOCITY  COMPONENTS  FOR  LORAN  C  (LC 

408)     48 

APPENDIX  C.     PROGRAM  DRIVLR  FORTRAN    51 

A.  GENERAL  REMARKS    51 

B.  DESCRIPTION  OF  PARAMETERS:      51 

APPENDIX  D.     TABLE  FOR  STANDARD  DEVIATIONS  AND 

COVARIANCES  OF  VELOCITY  COMPONENTS  FOR  MINI  RANGER  66 

APPENDIX  E.     AVERAGE  CURRENT  PROFILES    67 

APPENDIX  F.     CURRENT  PROFILES  FOR  SPECIFIC  TIMES    79 

LIST  OF  REFERENCES     85 

INITIAL  DISTRIBUTION  LIST      86 


LIST  OF  TABLES 

Table     1.  COMPARISON  OF  CLOCKS  WITH  THE  MINI  RANGER 

(MR)    2 

Table     2.  DIFFERENT  TIMES  BETWEEN  THE  CLOCKS  WITH  COM- 
MON TIME  BASE  (ADCP)     2 

Table     3.  LORAN  C  PHASE  CODES     13 

Table     4.  REFERENCE  STATION  SITES    30 

Table     5.  POSITION  ERROR  FOR  DIFFERENT  CROSS  ANGLES     .  .    32 
Table     6.  MAX-MIN  VALUES  OF  STANDARD  DEVIATIONS  AND 
COVARIANCES  FOR  VELOCITY  COMPONENTS  FOR  3 
MINUTE  TIME  INTERVALS  BETWEEN  POSITIONS  (MINI 

RANGER)     34 

Table     7.  SHIP  VELOCITY  COMPONENTS  FOR   MINI  RANGER. 

LORAN  (TD),  LORAN  (DISPLAY)  FROM  POINT  1  TO  POINT 

2  AT  3  MINUTE  TIME  INTERVAL 36 

Table     8.  SHIP  VELOCITY  COMPONENTS  FOR   MINI  RANGER, 

LORAN  (TD),  LORAN  (DISPLAY)  FROM  POINT  2  TO  POINT 
1  AT  3  MINUTE  TIME  INTERVAL     36 

Table     9.  COMPARISON  BETWEEN  MINI  RANGER  AND  LORAN 
(TD)  AND  THE  MINI  RANGER  AND  LORAN  (DISPLAY) 
FOR  LEGS  1  -  2  AND  2-1 37 

Table   10.  SHIP  VELOCITY  COMPONENTS  FOR   MINI  RANGER, 

LORAN  (TD).  LORAN  (DISPLAY)  FROM  POINT  1  TO  POINT 

3  AT  3  MINUTE  TIME  INTERVAL     38 

Table  11.  COMPARISON  BETWEEN  MINI  RANGER  AND  LORAN 

(TD)  AND  THE  MINI  RANGER  AND  LORAN  (DISPLAY) 

FOR  LEG  1-3 39 

Table   12.  SHIP  VELOCITY  COMPONENTS  FOR    MINI  RANGER. 
LORAN  (TD).  LORAN  (DISPLAY),  BOTTOM  TRACKING 


VI 


FROM  POINT  3  TO  POINT  2  AT  3  MINUTE  TIME  INTER- 
VAL         39 

Table  13.  SHIP  VELOCITY  COMPONENTS  FOR   MINI  RANGER, 
LORAN  (TD),  LORAN  (DISPLAY),  BOTTOM  TRACKING 
FROM  POINT  2  TO  POINT  3  AT  3  MINUTE  TIME  INTER- 
VAL         40 

Table  14.  COMPARISON  BETWEEN  MINI  RANGER  AND  LORAN 
(TD)  AND  THE  MINI  RANGER  AND  LORAN  (DISPLAY) 

FOR  LEGS  3  -  2  AND  2-3 40 

Table   15.  SHIP  VELOCITY  COMPONENTS  FOR   MINI  RANGER, 
LORAN  (TD).  LORAN  (DISPLAY),  BOTTOM  TRACKING 
FROM  POINT  3  TO  POINT  4  AT  3  MINUTE  TIME  INTER- 
VAL         41 

Table   16.  COMPARISON  BETWEEN  MINI  RANGER  AND  LORAN 
(TD)  AND  THE  MINI  RANGER  AND  LORAN  (DISPLAY) 

FOR  LEG  3-4 43 

Table  17.  AVERAGE  VELOCITY  COMPONENTS  FOR   MINI 

RANGER.  LORAN  (TD).  LORAN  (DISPLAY).  BOTTOM 

TRACKING     44 

Table  18.  TIME  INTERVAL  BETWEEN  POSITIONS  3  MINUTES  ...  48 
Table  19.  TIME  INTERVAL  BETWEEN  POSITIONS  6  MINUTES  ...  48 
Table  20.  TIME  INTERVAL  BETWEEN  POSITIONS  12  MINUTES  .  .  49 
Table  21.  TIME  INTERVAL  BETWEEN  POSITIONS  20  MINUTES  .  .  49 
Table  22.  TIME  INTERVAL  BETWEEN  POSITIONS  40  MINUTES     .  .    50 


Vll 


LIST  OF  FIGURES 

Figure     1.  CHART  OF  THE  WORK  AREA     3 

Figure     2.  VOLUME  SCATTERING  OF  SOUND    6 

Figure     3.  FIRST-SECOND  SOUND  DOPPLER  SHIFT    7 

Figure     4.  WATER  VELOCITY  VECTOR   8 

Figure     5.  VELOCITY  RESOLUTION  WITH  FOUR  ADCP  BEAMS    ...   8 

Figure     6.  SMALL  AND  LARGE  ERROR  VELOCITY     9 

Figure     7.  ADCP  AND  CURRENT  METERS    11 

Figure     8.  LOCATION  OF  WEST  COAST  LORAN  C  STATIONS    14 

Figure     9.  LORAN  C   GROUP  REPETITION  INTERVAL  (GRI)  AND 

TIME  DIFFERENCES  (TD)    15 

Figure   10.  THIRD  CYCLE  TRACKING  POINT    16 

Figure   11.  TYPICAL  HYPERBOLIC  FIX  GEOMETRY     17 

Figure   12.  HYPERBOLIC  GRADIENTS    19 

Figure   13.  HYPERBOLIC  GRADIENTS  AND  CROSSING  ANGLES     .    20 
Figure   14.  SPHERICAL  TRIANGLE  SPANAGEL-M-NORTH  POLE      .    21 

Figure  15.  PLANNING  THE  MOST  DESIRABLE  SITE    31 

Figure   16.  MAXIMUM  POSITION  ERROR     33 

Figure   17.  AVERAGE  CURRENT  PROFILE  FROM  POINT  1  TO 
POINT  3  USING  NAVIGATION  DATA  FROM  MINI 

RANGER     67 

Figure   18.  AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO 
POINT  2  USING  NAVIGATION  DATA  FROM  MINI 

RANGER     68 

Figure  19.  AVERAGE  CURRENT  PROFILE  FROM  POINT  2  TO 
POINT  3  USING  NAVIGATION  DATA  FROM  MINI 

RANGER     69 

Figure  20.  AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO 
POINT  4  USING  NAVIGATION  DATA  FROM  MINI 


VLU 


RANGER     70 

Figure  21.  AVERAGE  CURRENT  PROFILE  FROM  POINT  1  TO 

POINT  3  USING  NAVIGATION  DATA  FROM  LORAN  (TD)  71 

Figure  22.  AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO 

POINT  2  USING  NAVIGATION  DATA  FROM  LORAN  (TD)  72 

Figure  23.  AVERAGE  CURRENT  PROFILE  FROM  POINT  2  TO 

POINT  3  USING  NAVIGATION  DATA  FROM  LORAN  (TD)  73 

Figure  24.  AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO 

POINT  4  USING  NAVIGATION  DATA  FROM  LORAN  (TD)  74 

Figure  25.  AVERAGE  CURRENT  PROFILE  FROM  POINT  1  TO 
POINT  3  USING  NAVIGATION  DATA  FROM  LORAN 
(DISPLAY)     75 

Figure  26.  AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO 
POINT  2  USING  NAVIGATION  DATA  FROM  LORAN 
(DISPLAY)     76 

Figure  27.  AVERAGE  CURRENT  PROFILE  FROM  POINT  2  TO 
POINT  3  USING  NAVIGATION  DATA  FROM  LORAN 
(DISPLAY)     77 

Figure  28.  AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO 
POINT  4  USING  NAVIGATION  DATA  FROM  LORAN 
(DISPLAY)     78 

Figure  29.  CURRENT  PROFILE  FROM  POINT  1  TO  POINT  2  USING 
NAVIGATION  DATA  FROM  MINI  RANGER  FOR  A  SPE- 
CIFIC TIME  (14  33  1 1)    79 

Figure  30.  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  1  USING 
NAVIGATION  DATA  FROM  MINI  RANGER  FOR  A  SPE- 
CIFIC TIME  (15  06  10)    80 

Figure  31.  CURRENT  PROFILE  FROM  POINT  I  TO  POINT  2  USING 
NAVIGATION  DATA  FROM  LORAN  (TD)  FOR  A  SPE- 
CIFIC TIME  (14  33  1 1 )    81 

Fisure  32.  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  1  USING 


IX 


NAVIGATION  DATA  FROM  LORAN  (TD)  FOR  A  SPE- 
CIFIC TIME  (15  06  10)    82 

Figure  33.  CURRENT  PROFILE  FROM  POINT  1  TO  POINT  2  USING 
NAVIGATION  DATA  FROM  LORAN  (DISPLAY)  FOR  A 
SPECIFIC  TIME  (14  33  11)    83 

Figure  34.  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  1  USING 
NAVIGATION  DATA  FROM  LORAN  (DISPLAY)  FOR  A 
SPECIFIC  TIME  (15  06  10)    84 


H>'  , 


I.     INTRODUCTION 

The  relative  accuracy  of  navigation  data  is  very  important  in  the  processing 
of  Acoustic  Doppler  Current  Profiler  (ADCP)  data. 

The  shipboard  ADCP  is  an  instrument  which  measures  vertical  profiles  of 

■r 

currents.  The  raw  ADCP  data  is  collected  in  128  four  meter  bins,  and  is  used  in 
conjuction  with  a  navigation  system  from  which  the  position  of  the  ship  is  re- 
corded. 

In  this  experiment  we  will  try  to  evaluate  which  of  the  navigation  systems 
available  to  us  is  better  for  use  with  the  ADCP  in  order  to  determine  the  absolute 
current  velocity.  We  collected  data  with  MINI  RANGER,  LORAN  (TD), 
LORAN  (DISPLAY),  and  GPS  receivers.  Due  to  time  constraints,  however,  the 
GPS  data  has  not  been  included  in  the  analysis. 

Chapter  2  describes  the  cruise  from  which  the  data  analysed  in  this  thesis  is 
taken.  It  also  describes  the  time  synchonization  problems  which  existed  between 
the  various  types  of  navigational  equipment  used. 

Chapter  3  describes  the  ADCP  and  its  principles  of  operation. 

In  chapters  4  and  5  we  try  to  define  the  theoretical  random  error  for  the 
LORAN  (LC  40S)  and  the  MINT  RANGER  respectively. 

Chapter  6  gives  the  results  from  the  data  collected  in  a  cruise  of  Monterey 
Bay  while  Chapter  7  gives  a  summary  of  the  conclusions  from  this  study  together 
with  recommendations  for  future  work. 


II.     DATA  COLLECTION 

A.  CRUISE 

The  day  chosen  for  the  cruise  was  22  September  1989  with  the  ship  R/V  Point 
Sur.  During  the  cruise  the  following  equipment  belonging  to  the  NPGS  and  R/V 
Point  Sur  was  used: 

1.  Motorola  Mini  Ranger, 

2.  LORAN  C  receiver  (Internav  LC  408), 

3.  Trimble  GPS  receiver, 

4.  HP  Computer  logging  three  data  sources  simultaneously,  and 

5.  RDI  ADCP  with  the  Data  Acquisition  System  (DAS)  and  IBM  XT  com- 
patible computer. 

The  results  reported  in  later  sections  will  concern  only  the  MINI  RANGER 
and  LORAN  C  navigational  systems. 

The  ship  navigated  the  runs  from  point  1  to  point  2,  2  -  1,  1  -  3,  3  -  2,  2  -  3, 
and  3  -  4  (see  Figure  1).  The  speed  during  the  measurements  was  6  knots  which, 
from  previous  experience,  has  proven  effective  for  ADCP  data  collection.  In 
running  from  1  -  2  for  example,  a  turn  was  made  at  point  2  and  then  the  ship 
steamed  back  to  point  1  without  attempting  to  follow  the  identical  outward  track. 

The  average  depth  in  the  area  of  the  triangle  1  2  3  is  30  fathoms,  and  for  the 
run  3  -  4  it  varies  between  20  to  300  fathoms. 

B.  TIME  SYNCHRONIZATION 

During  the  cruise  the  four  pieces  of  equipment  used  each  had  clocks  which 
operated  independently  of  each  other.  In  order  to  synchronize  all  these  clocks,  at 
fixe  times  during  the  cruise  time  comparisons  were  made  between  them  by  using 
the  Mini  Ranger  as  the  common  time  base   (see  for  example  Table  1). 

These  times  were  then  converted  from  a  MR  time  base  to  an  ADCP  time 
base  (see  Table  2). 

The  standard  error  in  the  time  comparisons  was  at  least  1  second.  Once  a 
dedicated  multi-channel  data  logger  is  installed  which  allows  all  equipment  out- 


Figure  1.      CHART  OF  THE  WORK  AREA 


puts  to  be  recorded  on  a  common  time  base,  it  is  anticipated  that  this  synchroni- 
zation error  will  be  reduced  to  negligible  levels. 


Table  1.       COMPARISON  OF  CLOCKS  WITH  THE  MINI  RANGER  (MR) 


HP 

15:16:28 

MR 

15:16:28.6 

LC408 

15:14:00 

MR 

15:14:00.6 

ADCP 

15:15:10 

MR 

15:15:06.6 

Table  2.       DIFFERENT  TIMES  BETWEEN  THE  CLOCKS    WITH  COMMON 
TIME  BASE  (ADCP) 


HP 

15:16:28 

ADCP 

15:16:32.0 

LC40S 

15:14:00 

ADCP 

15:14:04.0 

MR 

15:15:06.6 

ADCP 

15:15:10 

With  five  tables  on  a  common  ADCP  time  base,  a  linear  regression  between 
the  different  times  was  performed  for  each  equipment  type  in  order  to  calculate 
the  coefficients  A  and  B  of  the  equation: 


Y=AX+ B 


(1) 


The  linear  clock  drifts  were  thus  calculated  and  used  to  assist  in  interpolating 
the  data  into  the  ADCP  time  base. 


III.     ACOUSTIC  DOPPLER  CURRENT  PROFILERS 

A.     INTRODUCTION 

The  shipboard  Acoustic  Doppler  Current  Profiler  (ADCP)  from  RD  Instru- 
ments is  a  state  of  the  art  instrument  which  measures  vertical  profiles  of  currents 
of  the  ocean  at  points  along  the  ship's  path  [RD  Instruments,  1989].  It  allows  the 
ocean  to  be  sampled  in  a  way  which  is  fundamentally  different  from  moored 
currents  or  drifters  (instruments  most  commonly  used  for  current  measurement). 

Consider  the  general  case  of  measuring  the  current  at  the  water  parcel  whose 
coordinates  are  X„(t)  using  an  instrument  whose  location  is  X0{t) 

XJt)  =  X0(t)  +  r  (2) 

Where  r  is  a  vector  between  the  two  locations  and 

r  =  Xw{t)  -  X0(t)  (3) 

From  (2)  and  (3) 

Xw(t)  =  X0(t)  +  [XM  -  X0(t)}  (4) 

The  current  U  at  the  measurement  point  XQ  +  r  will  be: 

r(Y        ,      dXw       dX0       d(Xw  -  Ap) 

^o  +  ')  =  — =  — + Jt (5) 

Equation  (5)  can  be  rewritten  as 

U(X0  +  r)  =  ^+  V{r)  (6) 

Where  V(r)  is  the  velocity  of  the  water  parcel  relative  to  the  instrument  posi- 
tion and  dXo.dt  is  the  velocity  of  the  instrument  itself  with  respect  to  the  Earth. 

Since  the  currents  are  a  difference  between  two  directly  measured  quantities 
dXo  dt  and  V.  the  measurement  of  currents  to  an  accuracy  of  1  cm  sec  from  a 


ship  travelling  at  10  knots  (approximately  5.2  m/sec  or  520  cm/sec)  requires  that 
both  dXo/dt  and  V  must  be  measured  to  an  accuracy  of  0.2%. 

B.     PRINCIPLES  OF  OPERATION  OF  THE  ACOUSTIC  DOPPLER 
CURRENT  PROFILER 

ADCP's  use  the  Doppler  effect  by  transmitting  a  succession  of  acoustic  pulses 
at  a  fixed  frequency  and  listening  to  the  resulting  backscattered  water  mass  ech- 
oes in  as  many  as  128  depth  cells  (bins)  over  a  depth  range  of  30  to  700  meters. 
When  scatterers  move  toward  the  ADCP,  the  echoes  heard  by  the  scatterers  is 
Doppler  shifted  to  a  higher  frequency.  The  amount  of  this  shift  is  proportional 
to  the  relative  velocity  between  the  ADCP  (ship)  and  scatterer  Figure  2. 
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Figure  2.      VOLUME  SCATTERING  OF  SOUND 


■6 


Part  of  these  Doppler-shifted  echoes  reflect  backwards  or  are  backscattered 
to  the  ADCP.  The  backscattered  echoes  appear  to  the  ADCP  as  if  the  scatterers 
were  the  source  of  the  echoes  (see  Figure  3). 

The  ADCP  hears  the  backscattered  echoes,  Doppler  shifted  a  second  time. 
Since  the  ADCP  both  transmits  and  receives,  the  Doppler  shift  is  doubled,  and 
the  equation  for  the  Doppler  shift  will  be 

FD  =  2Fs-^cosA  (7) 

where 
•     FD  is  the  Doppler  shift, 


MinH  Pulse    Moving  ScoUerers 


llllllr  ":f- 

Transducer 


Transducer       p^  gl    D0|>pler    Shift 

I 

ucer 

Second  Duppicr  Sliifl 


Transducer  ^TO, 


Figure  3.      FIRST-SECOND  SOUND  DOPPLER  SHIFT 

•  Fs  is  the  frequency  of  the  sound  when  everything  is  still, 

•  A  is  the  angle  between  the  acoustic  beam  and  the  water  velocity, 

•  V  is  the  relative  velocity  between  the  sound  source  and  the  sound  receiver, 
and 

.     C  is  the  speed  of  sound  (m/sec)  in  the  water  at  the  face  of  the  transducer, 
and  can  be  calculated  from  the  expression 

C  =  1449.2  +  4.67  -  5.51(T2r2  +  2.910'4r3       ^ 
+  (1.34  -  10-27%S  -  35)  +  1.61(T2£> 

where 

•  D  =  Depth,  in  meters, 

•  S  =  Salinity,  in  practical  salinity  units,  and 

•  T  =  Temperature,  in  degrees  Celsius. 

The  water  temperature  at  the  face  of  the  transducer  is  measured  by  the 
ADCP.  The  speed  of  sound  (C)  can  be  calculated  from  this  measured  water 
temperature  using  nominal  values  for  depth  (D)  and  salinity  (S). 
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WATER  VELOCITY  VECTOR 


Computer  analysis  of  the  Doppler  frequency  shift  of  backscattered  echoes 
from  each  bin  is  used  to  generate  a  precise  depth  segmented  picture  (or  profile) 
of  water  currents  throughout  the  water  column  bounded  by  the  path  of  the 
acoustic  beams. 

C.     THREE  DIMENSIONAL  CURRENT  VELOCITY  VECTORS 

The  ADCP  beams  each  measure  a  single  velocity  component,  i.e.  the  compo- 
nent of  velocity  toward  or  away  from  the  transduser.  When  using  multiple 
beams,  one  must  make  an  assumption  that  currents  are  the  same  (homogeneous) 
over  layers  of  constant  depth.  When  the  ADCP  uses  multiple  beams  pointed  in 
different  directions,  it  senses  different  velocity  components. 
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Figure  5.      VELOCITY  RESOLUTION  WITH  FOUR  ADCP  BEAMS 


Figure  5  shows  how  the  ADCP,  using  four  acoustic  beams,  computes  three 
velocity  components.  The  first  pair  of  beams  produces  one  horizontal  and  one 
vertical  component,  the  second  pair  produces  one  horizontal  (perpendicular  to  the 
first  horizontal)  and  again  one  vertical  component. 

The  product  SCALE  x  (Fj  +  V2  -  V3-  l\)  is  defined  as  the  error  velocity 
component,  where  Vx  to  F4  are  the  velocity  components  along  the  beam  directions 
Figure  6.  This  sum  should  be  close  to  zero.  The  error  velocity  component  allows 
us  to  evaluate  if  the  assumption  for  the  horizontal  homogenuity  is  reasonable. 
Also  this  is  another  way  to  estimate  data  quality. 
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Fisure  6.      SMALL  AND  LARGE  ERROR  VELOCITY 


In  Figure  6  there  are  two  different  situations,  in  the  first  the  velocity  is  the 
same  in  all  four  beams  at  a  constant  level,  in  the  second  the  velocity  in  one  beam 
is  different.  The  error  velocity  in  the  second  case  will  be  on  average  larger  than 
the  error  velocity  in  the  first  case.  The  difference  in  the  error  velocities  can  be 
from  different  currents  or  from  errors  caused  by  malfunctioning  equipment. 

D.     DEPLOYMENT  AREA 

The  RD  Instuments  ADCP  remotely  measures  water  flow  velocity  along  the 
lines  of  position  defined  by  the  four  narrow  vertically  inclined  acoustic  beams. 
To  insure  accurate  current  measurement  it  is  necessary  that  the  water  mass  in  the 


region  of  measurement  be  free  of  strong  acoustic  reflecting  objects  (e.g  platform 
members,  large  cables,  surface,  etc.)  within  a  ±15  degrees  conical  sector  along  the 
direction  of  each  of  the  beams.  However,  since  only  three  beams  are  required  for 
computation  of  three  axes  current  components,  in  applications  where  potential 
interference  objects  may  be  close  to  one  beam,  the  current  vector  may  be  calcu- 
lated from  the  other  three  beams. 

In  shallow  water  vertical  profiling  applications,  acoustic  interference  from  the 
surface  (upward  looking  ADCP)  or  bottom  (downward  looking  ADCP)  limit  the 
vertical  current  measurement  region  to  a  maximum  range  defined  by 

R{  max)  =  D  x  cos  <£  (9) 

where 

•  R(max)  =  maximum  profiling  range, 

•  D  =  distance  to  surface/bottom  boundary,  and 

•  </>  =  acoustic  beam  angle  relative  to  the  vertical. 

E.  BOTTOM  TRACKING  OPERATION 

The  ADCP  can  be  used  in  a  bottom  tracking  mode  to  give  direct  estimates 
of  the  velocity  components  of  the  ship.  In  the  experimental  data  collected  for  this 
thesis  we  used  the  bottom  tracking  on  for  the  runs  3  -  2,  2  -  3,  and  at  the  first  part 
of  the  run  3-4.  In  order  for  the  ADCP  to  gather  bottom  tracking  data,  bottom 
echoes  must  be  distinguished  from  other  echoes.  The  ADCP  transmits  a  dedi- 
cated bottom  track  ping  between  current  profiling  pings.  The  number  of  profiling 
pings  can  be  selected  by  the  user.  Consequently,  the  bottom  tracking  ping  data 
is  separated  from  the  current  profiling  data.  Bottom  echoes  are  identified  by  vir- 
tue of  their  greater  echo  strength  (the  bottom  returns  a  stronger  echo  than  the 
echo  from  a  profiling  bin). 

Ship  velocity  (relative  to  the  bottom)  is  measured  in  much  the  same  way  as 
current  velocity  (the  Doppler  shift  of  the  backscattered  bottom  echo  being  pro- 
portional to  the  ship  velocity).  Bottom  depth  is  detected  by  comparing  received 
echo  amplitude  to  a  detection  threshold  which  decreases  with  increasing  range  to 
the  bottom. 
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F.     VELOCITY  PROFILES 

The  most  important  feature  of  ADCP's  is  their  ability  to  measure  curreni 
profiles.  ADCP's  break  up  the  velocity  profile  into  uniform  segments  called  depth 
cells. 

Each  depth  cell  is  similar  to  a  single  current  meter,  but  there  are  two  basic 
differences.  The  first  difference  is  that  the  depth  cells  in  an  ADCP  are  always 
uniformly  spaced  while  current  meters  can  be  spaced  irregularly,  the  second  is 
that  the  ADCP  measures  average  velocity  over  the  depth  range  of  each  depth  cell 
while  the  current  meter  measures  current  only  at  the  current  meter  Figure  7. 
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Figure  7.      ADCP  AND  CURRENT  METERS 

The  advantage  of  the  ADCP  averaging  the  velocity  over  the  full  range  of  a 
depth  cell  is  that  it  reduces  the  effects  of  spatial  aliasing  (signals  at  frequencies 
higher  than  the  time  series  can  resolve  are  mistaken  for  low  frequency  signals). 

G.     ACOUSTIC  DOPPLER  CURRENT  PROFILER  DATA  PROCESSING 

The  raw  data  is  collected  in  128  four  meter  bins  over  a  three  minute  ensemble 
interval.  Generally  speaking,  either  a  MINI  RANGER,  LORAN  C  (calculated 
4>  ,  )■  from  time  differences  or  <f>  ,  X  recorded  from  the  display)  or  a  GPS  receiver 
are  used  for  navigation  and  a  position  recorded  every  30  seconds.  By  interpo- 
lation a  position  can  be  calculated  at  the  end  of  each  ensemble. 
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The  first  steps  in  processing  are  the  calculation  of  ship's  velocity  from  the 
navigation  data.  From  this  navigation  data  the  U  and  V  components  of  ship's 
velocity  are  calculated. 

The  next  step  in  processing  is  the  initial  determination  of  the  depth  to  which 
the  data  of  each  ensemble  remains  reliable.  For  our  case  the  7  ,  8  ,  9  bins  were 
chosen.  The  basic  criterion  comes  from  the  good  percent  of  return  echoes  (BIN 
STATISTICS  FILE).  By  subtracting  the  ship's  velocity  from  the  average  velocity 
within  the  chosen  reference  layer  an  absolute  reference  layer  velocity  for  each 
ensemble  is  obtained.  The  series  of  absolute  reference  velocities  is  then  filtered 
with  a  low  pass  Hamming  window  filter. 

Once  the  absolute  reference  velocity  is  determined  the  velocity  profiles  of  each 
ensemble  with  respect  to  the  reference  velocity  are  also  determined,  thus  yielding 
the  final  profiles  of  absolute  water  velocity.  The  remaining  profiles  of  absolute 
velocity  are  then  averaged  over  the  time  interval  for  each  run. 

In  appendix  E  there  are  18  average  current  profiles  from  different  navigation 
data. 
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IV.     THEORY  OF  LORAN   C 

A.  INTRODUCTION 

LORAN  C  (LONG  RANGE  NAVIGATION)  is  a  pulsed,  low  frequency 
(100  KHZ  carrier)  long-range  hyperbolic  navigation  system.  It  operates  on  the 
principle  that  the  difference  in  time  of  arrival  of  signals  from  two  stations,  ob- 
served at  a  point  in  the  coverage  area,  is  a  measure  of  the  difference  in  distance 
from  the  point  of  observation  to  each  of  the  stations. 

Measurements  of  the  ship's  latitude  and  longitude  (  4>  ,  X  )  and  Time  Differ- 
ences (TD)  from  a  LC  408  LORAN  C  receiver  were  recorded  throughout  the 
experiment  (22  September  1989)  at  intervals  of  30  seconds. 

B.  SYSTEM  DESCRIPTION 

LORAN  C  normally  requires  signals  from  at  least  a  Master  and  two  Sec- 
ondary stations  to  give  a  positional  fix. 

Transmitters  are  grouped  to  form  a  "chain"  of  which  one  station  is  labelled 
the  Master  (M)  and  the  others  are  called  secondaries  (X,Y,W)  as  shown  in  Figure 
8. 

The  chain  designators  for  LORAN  C  are  4  digit  numbers  which  indicate  the 
pulse  Group  Repetition  Interval  (GRI)  in  tens  of  micro-seconds.  For  example,  the 
west  coast  LORAN  C  chain  (Figure  8)  is  designated  9940  and  has  a  GRI  99400 
micro-seconds  (  see  Figure  9)  [LC  40S.  Operation  Manual] 

C.  PULSES-PHASE  AND  CODES-CYCLE  SELECTION 

Each  station  of  a  LORAN  C  chain  transmits  groups  of  pulses  (Table  3). 

Table  3.       LORAN  C  PHASE  CODES 


Master 

A  ++-+-+-  + 

B  4--+  +  +  +  +  - 

Secondary 

A  +  +  +  +  +-  + 

B  +  -+-++-. 
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Figure  8.      LOCATION  OF  WEST  COAST  LORAN  C  STATIONS 
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Figure  9.      LORAN  C    GROUP  REPETITION  INTERVAL  (GRI)  AND  TIME 
DIFFERENCES  (TD) 

The  Secondary  transmits  eight  pulses  and  the  Master  nine  (the  last  pulse  for 
indentification  and  blink  alarm).  The  pulses  are  phase  coded  to  improve  the  sig- 
nal to  noise  ratio  through  compression  and  to  distinguish  between  Master  and 
Secondary  transmissions. 

Low  frequency  radio  signals  propagate  over  the  earth's  surface  at  nearly  the 
velocity  of  light  in  a  vacuum.  However,  Maxwell's  equations  dictate  that  the 
ground  wave  (surface  wave)  will  be  influenced  slightly  by  the  surface  parameters 
of  geometry  and  electrical  properties. 

In  order  to  make  the  received  100  KHz  signal  more  stable  and  reliable  within 
a  given  coverage  area,  the  LORAN  C  radio  navigation  system  is  designed  as  a 
pulse  system  which  separates  the  ground  wave  from  the  skywave. 

Because  the  earth  parameters  remain  nearly  constant,  LORAN  C  has  dem- 
onstrated a  repeatability  of  quite  high  accuracy. 

The  high  accuracy  of  LORAN  C,  despite  long  ranges  from  transmitters,  is 
due  to  a  technique  called  "cycle  matching"  [Bowditch,  1984]. 

The  LC  408  tracks  the  third  cycle  cross-over  path  point  which  is  very  con- 
sistent between  transmissions  and  less  susceptible  to  skywave  interference  than 
later  cycles.  By  tracking  this  cycle  on  all  pulses,  high  accuracy  is  attained.  Figure 
10  shows  the  third  cycle  tracking  point  with  an  example  of  skywave  interference. 
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Figure  10.      THIRD  CYCLE  TRACKING  POINT 

D.     TIME  DIFFERENCE  EQUATIONS 

The  elapsed  time  between  the  arrival  of  the  signals  from  the  Master  and  a 
Secondary  stations  is  called  the  "Time  Difference"  (TD).  Each  observed  time 
difference,  or  rate,  provides  one  hyperbolic  line  of  position.  By  observing  the 
transmissions  from  four  stations,  three  hyperbolic  rates  are  measured  and  the 
position  can  be  determined  by  either  graphical  or  analytical  techniques. 

The  time  difference  observed  at  a  receiver  is  the  difference  in  arrival  times  for 
signals  from  the  master  and  one  secondary  transmitter  in  the  chain.  Because  all 
transmitters  share  the  same  frequencies,  their  signals  must  be  separated  in  time 
to  prevent  interference. 

The  chain  is  sychronized  so  that  the  Master  transmits  first  followed  by  each 
of  the  Secondaries.  The  transmission  of  each  secondary  is  specified  by  the  emis- 
sion delay  so  that  in  the  coverage  area  signals  from  one  station  will  overlap  an- 
other [Schenebele,  1979]. 

Suppose  that  at  a  point  P  there  is  a  receiver.  The  observed  time  difference 
using  the  Master  M  and  Secondaries  X,Y,W   are 
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Figure  11.      TYPICAL  HYPERBOLIC  FIX  GEOMETRY 


TDW=EDw+tw-tM 


TDK  =  EDX  +  tx  -  tj 


\f 


TDY  =  EDY+  tY-  r 


w 


(10) 

(11) 

(12) 


17 


where  EDW,  EDX,  EDY  are  the  emission  delay  for  W,  X,  Y.  The  terms  tw  ,tx,ty 
are  travel  times  from  the  W  ,  X  ,  Y  to  P.  The  term  tM  is  the  travel  time  from  the 
Master  M  to  P. 

In  order  to  express  the  time  difference  as  a  function  of  geographic  position  , 
the  travel  time  t  is  separated  into  additive  terms 

[  =  Jl-xD  +  f  (13) 

where: 

•  C  =  free  space  propagation  velocity, 

•  n  =  index  for  refraction  for  a  standard  atmosphere, 

•  D  =  geodetic  distance  from  the  transmitter  to  receiver,  and 

•  F  =  phase  factor  which  corrects  for  effects  of  the  earth's  surface  along  the 

path. 

Substituting  the  last  formula  into  the  equations  (10),  (11),  (12),  gives  the 
equations: 

TDW=EDw  +  ^x  {DW-DM)  +  FW-FM  (14) 

TDX  =  EDX  +  |x  (Dx  -  DM)  +  Fx  -  FM  (15) 

TDY  =  EDy  +  jr  x  (DY-DM)  +  FY-  FM  (16) 

All  these  equations  relate  the  time  differences  TDWJDX.TDY  to  the  dis- 
tances from  each  of  the  four  transmitters. 

Since  the  latitude  and  longitude  is  mathematically  calculated  from  the  time 
difference  numbers,  the  crossing  angles,  gradients,  and  signal  to  noise  ratios  are 
still  important. 

E.     HYPERBOLIC  GRADIENT,  CROSSING  ANGLES 

For  a  given  pair  Master,  Slave  a  family  of  hyperbolas  separated  by  a  costant 
value  can  be  plotted. 
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Figure  12.      HYPERBOLIC  GRADIENTS 

The  hyperbolas  in  Figures  12  and  13  are  all  separated  by  10  microseconds 
[Gann].  That  means  that  the  distance  on  the  baseline  represents  about  1852 
meters.  The  distance  between  successive  hyperbolas  increases  as  one  moves  to- 
wards the  baseline  extensions  or  away  from  the  baseline.  This  change  in  the  ac- 
curacy of  a  hyperbolic  Line  of  Position  (LOP)  that  occurs  relative  to  the  position 
of  the  LORAN  C  stations  is  known  as  the  gradient. 

The  angle  between  the  LOPs  is  called  the  crossing  angle  Figure  13.  The 
crossing  angle  must  be  bigger  than  30  degrees  and  smaller  than  150  degrees.  This 
fact  is  common  with  any  ranging  positioning  system,  but  is  further  limited  for  the 
LORAN  C  because  the  reception  of  three  stations  is  required  to  define  two  LOPs 
(direct  ranging  navigation  systems  provide  one  LOP  per  station). 
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Figure  13.      HYPERBOLIC  GRADIENTS  AND  CROSSING  ANGLES 

F.     ERROR  ANALYSIS  FOR  LORAN  C 

Any  discussion  of  accuracy  should  begin  with  the  definition  of  the  types  of 
accuracy. 

Repeatable  accuracy  (Repeatability,  Precision)  is  the  accuracy  with  which  a 
navigator  can  return  to  a  location,  the  coordinates  of  which  have  been  previously 
measured  using  the  same  system. 

Random  error  is  one  which  results  from  basic  limitations  in  the  method.  The 
characteristics  of  this  error  can  be  determined  by  statistical  analysis  of  a  sufficient 
number  of  measurements.  This  type  of  error  affects  the  repeatability.  These  er- 
rors are  identified  by  their  Gaussian  distribution. 

Absolute  accuracy  is  the  accuracy  with  which  a  navigator  can  determine  his 
point  position  in  terms  of  a  well  defined  coordinate  system  e.g.(</>  ,  X  ). 
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Systematic  errors  result  from  a  basic  (but  unrealized)  fault  in  the  method 
and  cause  the  values  to  be  consistently  biased  from  the  true  value.  It  affects  the 
absolute  accuracy,  and  often  cannot  be  detected  by  statistical  analysis. 

The  LC  408  microprocessor  uses  an  exacting  spherical  trigonometric  formula 
to  derive  its  latitude  and  longitude  positions.  Also  it  is  important  to  note  that  the 
LC  408  used  for  the  calculations  assumes  that  the  LORAN  signals  travel  over 
ideal  all  seawater  paths  from  the  transmitter  to  the  receiver.  Tp  correct  for  path 
over  land,  the  LC  408  allows  use  of  "additional  secondary  factors";  this  is  typi- 
cally used  only  in  small  regions  adjacent  to  known  fixed  locations. 

For  the  error  analysis  given  shortly,  the  shape  of  the  earth  is  approximated 
by  a  sphere,  and  it  is  assumed  that  all  the  west  longitudes  are  positive  . 

In  order  to  find  the  standard  deviations  in  the  velocities,  we  established  a 
fixed  point  in  the  roof  of  Spanagel  Hall  to  enable  an  assessement  of  the  standard 
deviations  in  the  three  time  differences  and  after  it  to  propagate  the  error  at  the 
positions  and  the  velocities. 

Consider  the  spherical  triangle  PSM  (see  Figure  14  [Cross,  1981]). 


Figure  14.      SPHERICAL  TRIANGLE  SPANAGEL-M-NORTH  POLE 
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Where  (M)  can  be  the  Master,  S  can  be  the  point  on  the  Spanagel  roof,  and 
P  can  be  the  North  pole.  We  try  to  find  an  expression  for  the  distance  between 
Master  and  Spanagel  roof  with  respect  to  <j>  ,  /  (coordinates  of  the  point  on  the 
Spanagel  roof), 

m  =  90 -<j)  (17) 

s  =  90-4>M  (18) 

P  =  X-?.M  (19) 

where  <f>M  and  XM  are  the  coordinates  of  the  Master  station.    Using  the  law  of 
cosines  for  the  side  p, 

cos  p  =  cos  m  cos  s  +  sin  m  sin  s  cos  P  (20) 

Substitute  (17)  ,  (18)  ,  (19)  into  (20)  to  yield 

cos  p  =  sin  4>M  sin  <f>  +  cos  <f>M  cos  4>  cos(>.  -  AM)  (2 1 ) 

Differentiate  each  part  of  (21) 

—  sin  pdp  =  cos  4>m  sin  4>d<j)^  +  sin  4>m  cos  4>d4> 

—  sin  <f>  u  cos  <f>  cos(/  —  }-\t)d<t>M 

—  cos  (f>M  sin  <j)  cos(/  —  /.^)d4>  (22) 

—  cos  4>Af  cos  (/>  sin(/  —  XM)dA 
+  cos  </)iU  cos  4>  sin(/  —  /W)tf/A/ 

Since  the  position  of  Master  station  is  a  fixed  point  d).M  =  d(f>M  =  0  equation  (22) 
becomes 

—  sin  pdp  =  [  sin  <f>M  cos  <f>  —  cos  <f>  u  sin  <^>  cos(/  —  /iV/)W 

-  cos  </)M  cos  (f>  sin(/  —  /  v/)t// 

We  now  introduce  into  each  term  in  the  right  part  of  (23)  the  factor  sin  p.   Using 
the  relationship  between  two  angles  and  three  sides 
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sin  p  cos  5  =  cos  s  sin  m  -  sin  s  cos  m  cos  P  (24) 

Substituting  for  s,  m  and  P  as  appropriate,  we  have 

sin  4>m  cos  4>  ~~  cos  4>m  sin  <f>  cos(/  —  Aw)  = 
cos(90  -  <f>M)  sin(90  -</>)-  sin(90  -  4>M)  cos(90  -  </>)  cos(/  -  AM)  =      (25) 

sin/>  cos  5 

From  the  law  of  sines  we  know  that 

sin  m        sin  p 


(26) 


sin  .1/       sin  P 
OR 

sin  m  sin  P  =  sin  p  sin  A/  (28) 

Substituting  for  m  and  P  and  multiplying  both  sides  by  cos  <f>M  ,  we  have 

cos  4> .;/  cos  (f)  sin(/  —  /.  u)  = 

cos  4>X!  sin(90  -  0)  sin(/  -  AM)  =  (29) 

sin  p  sin  M  cos  <f>M 

Substituting  (25)  and  (29)  into  (23).  we  have 

dp  =  —  cos  Sd(p  +  sin  M  cos  (fr^d/.  (30) 

Actually  in  order  to  convert  the  angle  dp  into  distance  da  (arc  length),  using 
spherical  approximation  multiply  by  the  mean  radius  of  the  earth  R 

da  =  -  R  cos  Sd<j)  +  R  sin  M  cos  <p\fd/.  (31) 

From  the  last  equation  we  have 

=  ->RcosS  (32) 


C(p 


23 


da 
-—-  =  R  cos  4>M  sin  M  (33) 

where  M  and  S  are  the  interior  angles  of  the  triangle  PSM.    In  similar  way  can 
form  the  equations  for  the  other  spherical  triangles 

— -  =  -Rcos  SY  (34) 

ocp 

— —  =  -  R  cos  4>x  sin  X  (35) 

c  /. 

The  sign  changes  in  the  second  equation  due  to  different  geometry 

da-) 


c<t> 
ca2 


=  -  R  cos  S2  (36) 


CA 


=  Rcos<f)Ysm  Y  (37) 


C<7: 

tt  =  -  R  cos  5,  (38) 

00 


-^r-=/?cos0™sin  ir  (39) 

c/ 


Next  let 

4  =  *((M)  -  *i(<M)  (40) 

rf2  =  *(<M)  -  *2(<M)  (41) 

rf3  =  *(</>,  ;.)  -  «3(0,  /)  (42) 

Suppose  that  the  velocity  of  the  the  electromagnetic  propagation  is  constant  and 
equal  (C)  then: 
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\  = 

4 

C 

\= 

d2 
C 

V= 

d3 
C 

(43) 


(44) 


(45) 


We  propagate  the  errors  into  these  expressions  by  forming  the  Jacobian  ma- 
trix. We  assume  here  that  time  is  subject  to  negligible  error.   This  matrix  is  given 

by 


./ 


c 


dcf) 

Cc12 

c4> 


dd 


CA 

cd3 
a  A 


(46) 


and 


dd{         cd{(i),A)        cax{4>. /.) 

d<p  c</>  c4> 

=  -  R  cos  S  +  R  cos  Si 

=  R(  cos  5j  -  cos  5) 
da{<j>.  /)        da2(<t>,  a) 


cd2 

c4> 


dd3 

C(j) 


c(b  ccf) 

-  R  cos  S  +  R  cos  S2 

R(  cos  S2  —  cos  S) 
ca(<f).  a)       da3(<j>,  A) 


c<i>  c<f) 

=  —  R  cos  S  +  R  cos  S3 

=  R(  cos  S3  -  cos  S) 


(47) 


(48) 


(49) 


25 


ddx         da{4>,?.)        dax{4>,)) 

dA  cX  dk 

=  R  cos  4>M  sin  M  +  R  cos  <j>x  sin  X 

=  R{  cos  <pM  sin  M  +  cos  <f>x  sin  JQ 

dd2         da(<t), ).)        da2{4>,  /) 

3  i 

<£rsin  F 

M  —  cos  0  F  < 


3/  3A  3/ 

=  ^  cos  <£M  sin  A/  —  R  cos  </>F  sin  Y 

=  #(  cos  <£jVf  sin  M  —  cos  0Ksin  Y) 
dd3         da{4), ).) 

os  <f>M  sin  M  —  R  cos  </>^sin  IV 


VA  US. 

=  /?  cos  4>M  sin  A/  —  R  cos  $  ^  sin  IV 
=  R{  cos  </> M  sin  A/  -  cos  <j>wsm  W) 


(50) 


(51) 


(52) 


C  =  3  x  108(A//  sec)  =  3  x  10-1(A///zjec) 


Now  define  the  variance  covariance  matrix  of  LORAN  time  differences,  as- 
suming that  there  is  no  correlation  between  the  three  time  differences. 


V  _ 


0 

2 

0 


0 


(54) 


The  LC  408  is  provided  with  control  of  the  filtering  time  constant.  Selectable 
time  constants  are  0  seconds  for  no  filtering.  5  seconds  for  good  response,  10 
seconds  for  slower  response,  and  20  seconds  for  slow  craft  or  monitor  applica- 
tions. 

In  order  to  measure  the  standard  deviations  for  the  time  differences,  we  re- 
corded data  for  four  days  on  Spanagel  Roof,  each  day  (24  hours)  with  a  different 
time  constant.  This  resulted  in  four  different  variance  covariance  matrices. 
There  are  all  given  in  Appendix  A. 

From  these  results  we  conclude  that  the  filters  don't  work  for  the  time  dif- 
ference between  Master  and  Secondary  Whisky  because  the  standard  deviation 
remains  approximately  the  same  for  the  different  time  filters.  We  have  no  ex- 
planation for  this  beyond  it  being  an  instrumental  problem.    We  also  conclude 
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that  for  monitoring  applications,  the  longer  the  filtering  time,  the  smaller  the 
standard  deviations. 

The  weight  matrix  P  =  Ir1  is  given  by 


P  = 


Mai 

0 

0 


0 

Wo) 
0 


0 
0 


(55) 


The  variance  covariance  matrix  for  the  position  [Uotila,  1986]  will  be: 


L(4>,  A)  =  [JT  x  P  x  J) 


-l 


(56) 


Define  the  matrix  In  as 


-c 


,0  0 


(57) 


Now  we  form  the  block  matrix  1  (  </>]/], ^2/2)  which  is  the  variance 
covariance  matrix  for  two  positions  (the  error  characterises  are  the  same  for  the 
two  positions).  If  it  is  assumed  that  there  is  no  correlation  between  the  two  posi- 
tions, then 


£(<Ml,4>2/2)  = 


^l./l) 


v 
"0 


v 
^0 


!(02>  'J 


(58) 


In  order  to  convert  the  difference  in  latitude  and  difference  in  longitude  (radians) 
between  two  positions  in  meters,  apply 


A  =  A<f>  x  R 


(59) 


B  =  A/  x  cos  <j)m  x  R 


(60) 


where  <j>m  is  the  mean  latitude  of  the  work  area,  and   R  is  the  mean  radius  oi  the 
earth  R  =  6371  KM.    Let 


C,  =/? 


(61) 
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C2  =  cos  4>m*  R 


(62) 


The  components  of  the  ship  velocity  V  ,  U  (Northing,  Easting)  are  then  given  by: 


At        At      KH2      ^v 


(63) 


U  = 


B   _Ci 
At       At 


x  (;.2  - ;.}) 


(64) 


The  Jacobian  matrix  for  the  velocity  will  be 


/  = 


At 


di 


d<t>i 
cl1 

c4i 


zi/ 


d\ 


eu 


d/. 


dV 

d<p2 

cU 

C$2 


3  17 

cv 
dX2 

dU 


(65) 


and 


/  = 


1 

-Ci 

0 

Ci 

0 

At 

0 

-Q 

0 

c2 

(66) 


Consider  time  intervals  between  positions   3  ,  6  .  12  ,  20  and  40  minutes.  The 
variance  covariance  matrix  for  the  velocity  components  will  be: 


'V,  u 


(67) 


Appendix  B  gives  5  tables  for  standard  deviations  and  covariances  of  velocity 
components  with  different  time  intervals  between  positions.  From  these  tables 
we  see  clearly  that  as  the  time  interval  between  the  positions  increases  the  stand- 
ard deviations  for  the  V  and  U  components  decreases  (assuming  constant  course 
and  speed  of  the  ship  for  the  given  interval).  Also  we  note  that  when  time  aver- 
aging over  20  minutes  or  longer,  the  choice  of  filter  is  not  critical. 

Appendix  C  gives  a  program  (DRIVLR  FORTRAN)  written  by  the  NOAA 
which  calculates  the  geographic  position  of  the  ship  from  two  LORAN  Time 
Differences  (TD)  and  a  Dead  Reckoning  (DR)  position  of  the  ship  (inverse  com- 
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putation).  It  can  also  be  used  to  convert  the  geographic  position  of  the  ship  into 
LORAN  Time  Differences  (forward  computation). 
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V.     THEORY  OF  MINI  RANGER 

A.  INTRODUCTION 

The  Mini  Ranger  Falcon  was  used  during  the  cruise  for  navigational  posi- 
tioning. The  standard  Mini  Ranger  operates  at  up  to  37  kilometers  (about  20 
nautical  miles)  with  a  probable  range  measurement  error  of  2  meters. 

It  operates  at  microwave  frequencies  and  requires  that  line  of  sight  be  main- 
tained between  the  reference  stations  and  the  receiver  transmitter.  Significant 
obstuctions  such  as  land  masses,  buildings,  or  dense  foliage  will  interfere  the  op- 
eration of  the  system. 

B.  REFERENCE  STATIONS 

The  positions  of  the  reference  station  sites  in  the  UTM  coordinate  system  are 
[Krioneritis.  1989]  listed  in  Table  4. 

Table  4.       REFERENCE  STATION  SITES 


STATION 

X  EASTING 

Y  NORTHING 

CODE 

TREVOR 

585260.161  M 

4092490.284  M 

15 

PACK 

609863.128  M 

4076611.345  M 

1 

HAYES 

607621.289  M 

4055915.264  M 

12 

C.     RANGE  POSITIONING 

The  Mini  Ranger  determines  a  two  range  position  when  the  lengths  of  all 
three  sides  of  a  triangle  are  known  (trilateration).  Suppose  that  a  certain  task  re- 
quires activity  between  points  C  .  D  ,  E  (Figure  15),  the  reference  stations  are  at 
points  A  and  B.  in  the  acceptable  area  the  angle  of  intersection  between  the  two 
rancc  lines  is  between  30  and  150  degrees. 
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Figure  15.      PLANNING  THE  MOST  DESIRABLE  SITE 

Arcs  AXB  and  AYB  define  the  maximum  and  minimum  distance  to  maintain 
between  30  and  150  degrees  respectively. 

Generally  the  Mini  Ranger  uses  the  method  of  least  squares  to  calculate  po- 
sitions from  three  or  more  reference  stations.  Operation  in  areas  where  the  ge- 
ometry is  poor  will  degrade  positional  accuracy  [Mini  Ranger  operation  manual, 
1981]. 

D.     ERROR  ANALYSIS  FOR  MINI  RANGER 

While  the  Mini  Ranger  measures  ranges  well,  the  actual  position  error  re- 
sulting from  a  set  of  range  readings  depends  upon  a  number  of  contributing  fac- 
tors: 

1 .  system  error,  and 

2.  geometry  changes. 
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Figure  16  shows  the  maximum  positional  error  for  two  range  geometries  of  30  , 
90  and  150  degrees  [Mini  Ranger  operation  manual,  1981]. 


POSITION  ERROR  = 


STANDARD  ERROR 

sin  6 


(68) 


The  position  error  for  various  two  range  geometries,  assuming  range  error  of  a 
single  measurement  of  2  meters  are  given  in  Table  5. 


Table  5.       POSITION    ERROR    FOR    DIFFERENT 
CROSS  ANGLES 


CROSSING  ANGLE 

POSITION  ERROR 

150  degrees 

7.7  meters 

1 20  degrees 

4.0  meters 

90  degrees 

2.8  meters 

60  degrees 

4.0  meters 

30  degrees 

7.7  meters 

For  this  experiment  we'll  take  the  information  for  the  variance  covariance 
matrix  of  the  positions  from  the  Mini  Ranger  data  processing  (see,  for  example, 
Krioneritis,  19S9).    Define  the  matrix  Z0  as 


v     ro       o] 

^  ~  Lo        oj 


(69) 


The  variance  covariance  matrix  for  the  position  will  be 


2W)  = 


tx,y 


ffX,Y 

2 


(70) 


We  now  form  the  block  matrix  I  (  Xu YuX2t Y2)  which  is  the  variance 
covariance  matrix  for  two  positions.  If  it  is  assumed  there  is  no  correlation  be- 
tween the  two  positions,  then 
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Figure  16.      MAXIMUM  POSITION  ERROR 


Z(x1y1x2y2)  = 


s(A'„  rj 


i(A2,  r2) 


71) 


Define  the  velocity  components  as 


V  = 


y2  -  *i 

Ar 


72) 


t/  = 


A  2  ~  Xi 

At 


(73) 
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As  we  did  with  LORAN  C,  we  propagate  errors  into  these  expressions  by  forming 
the  jacobian  matrix  for  the  velocity  components.  We  again  assume  that  time  is 
subject  to  negligible  error.   This  matrix  is  given  by 


/  = 


At 


dV 
dX} 
dU 
dX-i 


dV 

ILL 


dV 

dx2 

ML 
ex. 


21/ 


dl 


dY2 

ILL 

dYi 


(74) 


OR 


J  = 


At 


0 
-1 


-1 

0 


(76) 


The  variance  covariance  matrix  for  the  velocity  components  will  be 


^f,  u  -  J  x  ^x1y1x2y2  x  J' 


(77) 


Table  6.  MAX-MIN  VALUES  OF  STANDARD  DEVIATIONS  AND 
COVARIANCES  FOR  VELOCITY  COMPONENTS  FOR  3  MINUTE 
TIME  INTERVALS  BETWEEN  POSITIONS  (MINI  RANGER) 


MAX-MIN 
VALUES 

Oy 

<*u 

av,u 

MAXIMUM 

1.6  cm,  sec 

1.8  cm,  sec 

2.Scm2/  sec2 

MINIMUM 

1 .2  cm  sec 

1.3  cm  sec 

0.2cm2/  see2 

Table  (6)  shows  that  the  the  standard  error  in  the  velocity  components 
doesn't  vary  greatly  given  the  geometry  in  the  work  area.  The  results  come  from 
6  different  points  in  the  whole  work  area. 

Appendix  D  gives  the  velocity  error  results  for  different  time  intervals.  Simi- 
lar to  the  LORAN  (Appendix  B),  the  standard  deviations  for  the  velocity  com- 
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ponents    decrease    when    the    time    interval    between    the    positions    increases 
(assuming  constant  course  and  speed  of  the  ship  for  the  given  time  interval). 

Comparing  the  theoretical  results  from  Appendices  B  and  D,  we  can  see  that 
the  MINI  RANGER  is  not  only  a  more  stable  navigation  system  than  the 
LORAN  C,  but  also  that  it  provides  in  3  minutes  the  level  of  accuracy  in  the  U 
and  V  velocity  components  w7hich  require  20  minutes  to  obtain  with  LORAN  C. 
Because  of  the  practical  difficulties  in  maintaining  a  constant  ship's  heading  and 
speed  for  20  minutes,  and  thus  the  need  for  an  averaging  process  over  such  a  time 
interval,  the  value  of  using  a  MINI  RANGER  wherever  possible  is  even  further 
enhanced. 
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VI.     RESULTS 

In  this  chapter  there  are  eleven  Tables  which  summarise  the  results  of  the 
computation  of  the  ship  velocity  components  from  MINI  RANGER,  LORAN 
Time  Differences  (TD),  LORAN  (DISPLAY)  and  Bottom  Tracking.  We  choose 
to  separate  the  LORAN  (TD),  and  LORAN  (DISPLAY)  results  because  of  the 
rounding  off  which  occurs  when  using  the  LORAN  (DISPLAY). 

The  LORAN  results  given  here  were  obtained  using  different  filters  (2  -  1,  1 
-  3  filter  with  time  constant  5  seconds,  1  -  2  filter  with  time  constant  20  seconds, 
3  -  2,  2  -  3  filter  with  time  constant  10  seconds,  3  -  4  filter  with  time  constant  0 
seconds  ). 

Table  7.  SHIP  VELOCITY  COMPONENTS  FOR  MINI  RANGER.  LORAN 
(TD),  LORAN  (DISPLAY)  FROM  POINT  1  TO  POINT  2  AT  3  MIN- 
UTE TIME  INTERVAL. 


TIME 

MINI  RANGER 

LORAN  (TD) 

LORAN  (DIS- 
PLAY) 

U 
(cm  sec) 

V 

(cm  sec) 

U 
(cm  sec) 

V 

(cm  sec) 

U 
(cm  sec) 

V 

(cm  sec) 

14  24  11 

-126.5 

-277.1 

-125.4 

-273.5 

-126.2 

-263.3 

14  27  10 

-125.9 

-276.0 

-123.2 

-274.8 

-124.0 

-286.4 

14  30  10 

-123.3 

-277.6 

-120.6 

-278.1 

-115.7 

-271.5 

14  33  11 

-124.4 

-280.5 

-124.4 

-274.6 

-124.0 

-282.8 

14  36  11 

-106.3 

-282.8 

-107.4 

-282.5 

-109.3 

-281.6 

14  39  11 

-108.7 

-285.1 

-103.0 

-296.2 

-105.2 

-275.6 

14  42  10 

-109.3 

-280.7 

-116.5 

-284.5 

-113.4 

-268.4 
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Table  8.  SHIP  VELOCITY  COMPONENTS  FOR  MINI  RANGER,  LORAN 
(TD),  LORAN  (DISPLAY)  FROM  POINT  2  TO  POINT  1  AT  3  MIN- 
UTE TIME  INTERVAL 


TIME 

MINI  RANGER 

LORAN  (TD) 

LORAN  (DIS- 
PLAY) 

U 
(cm  sec) 

V 
(cm /sec) 

U 
(cm/sec) 

V 
(cm/sec) 

U 
(cm/sec) 

V 

(cm,  sec) 

14  57  10 

135.0 

265.3 

134.3 

254.5 

135.0 

258.0 

15  00  10 

133.7 

261.3 

134.4 

256.1 

132.2 

240.7 

15  03  10 

13S.9 

255.4 

134.7 

253.7 

138.1 

265.0 

15  06  10 

154.4 

252.0 

154.0 

245.6 

148.7 

244.1 

15  09  10 

155.9 

253.1 

154.7 

256.5 

151.7 

247.9 

15  12  10 

148.8 

255.4 

146.2 

252.3 

154.1 

261.3 

15  15  10 

155.7 

251.8 

154.9 

250.7 

143.2 

237.2 

From  Tables  7  and  8  we  can  see  a  discontinuity  beginning  at  14  36  11  near 
the  end  of  leg  1  -  2  and  ending  at  15  03  10  at  the  beginning  of  leg  2  -  1.  This 
discontinuity  occurs  at  the  same  position  down  each  leg  and  is  of  approximately 
the  same  magnitude  in  each  case.  While  it  is  tempting  to  view  this  discontinuity 
as  being  due  to  a  surface  current,  the  inconsistency  in  sign  prohibits  such  an  ex- 
planation. It  can  only  be  concluded  therefore  that  it  is  due  to  a  ship  navigation 
correction. 

In  these  tables  we  see  clearly  that  the  LORAN  (DISPLAY)  results  show  a 
much  geater  variability  than  the  LORAN  (TD).  Also  if  we  calculate  the  mean 
value  and  the  standard  deviation  of  the  mean  of  the  absolute  difference  of  the  U 
and  V  components  between  the  MINI  RANGER  and  LORAN  (TD).  and  the 
MINI  RANGER  and  LORAN  (DISPLAY),  we  have  Table  9. 
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Table  9.  COMPARISON  BETWEEN  MINI  RANGER  AND  LORAN  (TD)  AND 
THE  MINI  RANGER  AND  LORAN  (DISPLAY)  FOR  LEGS  1  -  2  AND 

2-  1. 


LEG 

MINI  RANGER  -  LORAN  (TD) 

MINI  RANGER  -LORAN  (DISPLAY) 

A  U  (frAV) 

A  V  (<rAK) 

A  U  (<rAl7) 

A  V  (<rAK) 

1  -2 

2.9(0.9) 

4.5(2.1) 

2.9(0.9) 

7.9(1.8) 

2  -  1 

1.5(0.5) 

4.5(1.2) 

4.2(1.6) 

10.2(2.1) 

After  using  the  FISHER-BEHERENS  test  [Hamilton,  1964]  to  compare  Ac/  from 
MINI  RANGER-LORAN  (TD)  against  At/  from  MINI  RANGER-LORAN 
(DISPLAY)  and  AV  against  AV  we  see  that  these  are  statistically  different  at 
95%  confidence  interval  in  3  out  of  4  cases.  This  provides  strong  evidence  to 
suggest  that  the  results  from  LORAN  (TD)  and  LORAN  (DISPLAY)  are  from 
different  statstical  populations. 

Table  10.  SHIP  VELOCITY  COMPONENTS  FOR  MINI  RANGER,  LORAN 
(TD),  LORAN  (DISPLAY)  FROM  POINT  1  TO  POINT  3  AT  3 
MINUTE  TIME  INTERVAL 


TIME 

MINI  RANGER 

LORAN  (TD) 

LORAN  (DIS- 
PLAY) 

U 
(cm  sec) 

V 
(cm  sec) 

U 
(cm/sec) 

V 

(cm  sec) 

U 
(cm/sec) 

V 
(cm  sec) 

15  30  10 

-272.4 

-116.1 

-272.0 

-118.6 

-272.8 

-123.8 

15  33  11 

-273.1 

-115.8 

-274.6 

-116.5 

-253.5 

-106.8 

15  36  11 

-277.4 

-117.8 

-276.6 

-116.1 

-291.8 

-126.8 

15  39  10 

-277.1 

-117.9 

-276.8 

-121.4 

-273.7 

-104.1 

15  42  10 

-268.4 

-122.7 

-269.5 

-123.2 

-274.8 

-130.6 

15  45  10 

-266.7 

-126.4 

-268.6 

-130.2 

-264.4 

-134.0 
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15  48  10 

-265.9 

-123.9 

-264.7 

-119.4 

-264.7 

-123.9 

15  51  10 

-268.9 

-126.7 

-265.9 

-119.8 

-264.3 

-113.5 

15  54  10 

-269.8 

-123.9 

-263.8 

-118.7 

-276.3 

-123.9 

15  57  10 

-266.7 

-126.1 

-267.2 

-124.8 

-250.8 

-125.8 

16  00  10 

-269.5 

-120.5 

-269.5 

-114.8 

-269.9 

-101.3 

16  03  11 

-265.9 

-127.2 

-259.5 

-127.9 

-273.8 

-140.1 

16  06  10 

-268.6 

-125.5 

-267.4 

-125.2 

-260.6 

-123.8 

Table  11.  COMPARISON  BETWEEN  MINI  RANGER  AND  LORAN  (TD)  AND 
THE  MINI  RANGER  AND  LORAN  (DISPLAY)  FOR  LEG  1  -  3. 


LEG 

MINI  RANGER  -  LORAN  (TD) 

MINI  RANGER  -LORAN  (DISPLAY) 

A  L  (<tm ■) 

A  V  (<rA„) 

A  L   (<7A£/) 

A  V  (au) 

1  -  3 

1.9(0.6) 

2.9(0.6) 

7.0(1.7) 

7.9(1.-) 

Tabic  10  shows  similar  results  for  leg  1  -  3.  From  Table  11  it  is  also  obvious 
that  the  results  which  comes  from  the  LORAN  (TD)  and  the  LORAN  (DIS- 
PLAY) are  different  for  leg  1  -  3. 

Table  12.  SHIP  YELOCITY  COMPONENTS  FOR  MINI  RANGER.  LORAN 
(TD),  LORAN  (DISPLAY).  BOTTOM  TRACKING  FROM  POINT  3  TO 
POINT  2  AT  3  MINUTE  TIME  INTERVAL 


TIME 

MINI  RANGER 

LORAN  (I'D) 

LORAN  (DIS- 
PLAY) 

BOTTOM 
TRACKING 

L 

(cm  sec) 

V 

(cm  sec) 

L 

(cm  sec) 

V 

(cm  sec) 

L 

(cm  sec) 

V 

(cm  sec) 

L 

(cm  sec) 

V 

(cm  sec) 

16  16  18 

272.4 

-103.5 

267.9 

-104.4 

259.7 

-105.4 

263.5 

-116.5 

16  19  18 

277.9 

-93.2 

280.2 

-94.6 

274.1 

-84.0 

271.4 

-103.5 

16  22  18 

275.5 

-95.3 

272.1 

-95.5 

279.7 

-111.8 

268.5 

-106.0 

16  25  IS 

2"". 3 

-99.2 

275.1 

-97.5 

276.6 

-SO. 7 

571     3 

-110.1 
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16  28  17 

278.0 

-98.2 

282.2 

-95.4 

275.7 

-98.7 

272.4 

-108.1 

16  31  19 

280.1 

-106.6 

279.8 

-97.7 

281.3 

-99.1 

274.2 

-109.5 

Table  13.  SHIP  VELOCITY  COMPONENTS  FOR  MINI  RANGER,  LORAN 
(TD),  LORAN  (DISPLAY),  BOTTOM  TRACKING  FROM  POINT  2  TO 
POINT  3  AT  3  MINUTE  TIME  INTERVAL 


TIME 

MINI  RANGER 

LORAN  (TD) 

LORAN  (DIS- 
PLAY) 

BOTTOM 
TRACKING 

L 

(cm  sec) 

V 

(cm  sec) 

U 

(cm  sec) 

V 

(cm  sec) 

U 

(cm  sec) 

V 

(cm  sec) 

L 

(cm  sec) 

V 

(cm  sec) 

16  52  IS 

-260.3 

124.3 

-256.1 

124.3 

-256.5 

124.3 

-253.0 

128.5 

16  55  17 

-259.2 

118.4 

-261.1 

114.4 

-271.0 

121.5 

-252.8 

124.8 

16  58  17 

-262.8 

116.4 

-264.7 

115.9 

-249.9 

108.1 

-257.1 

125.9 

17  01  17 

-269.5 

122.7 

-266.1 

118.0 

-277.1 

126.5 

-262.8 

12S.4 

17  04  18 

-271.1 

117.8 

-268.9 

120.9 

-266.6 

118.5 

-264.9 

125.1 

17  07  18 

-264.2 

116.4 

-265.3 

108.3 

-266.8 

114.2 

-258. 8 

123.0 

17  10  IS 

-267.0 

1 1 5.4 

-263.6 

119.4 

-260.2 

120.8 

-262.9 

123.2 

17  13  17 

-264.0 

121.9 

-260.9 

120.7 

-264.3 

116.0 

-258.4 

128.0 

From  Tables  12  and  13  we  can  see  that  the  more  consistent  results  comes 
from  MINI  RANGER.  The  variability  of  the  LORAN  (TD)  results  are  smaller 
than  the  variability  of  the  LORAN  (DISPLAY)  results.  We  also  see  that  the 
Bottom  Tracking,  while  giving  values  to  the  same  internal  consistency  as  the 
MINI  RANGER,  gives  results  which  are  systematically  biased. 

Table  14.  COMPARISON  BETWEEN  MINI  RANGER  AND  LORAN  (TD)  AND 
THE  MINI  RANGER  AND  LORAN  (DISPLAY!  FOR  LEGS  3  -  2  AND 


LEG 

MINI  RANGER-  LORAN 

(TD) 

MINI  RANGER  -LORAN 
(DISPLAY) 

MINI  RANGER  -  BOT- 
TOM TRACKING 

A  L   (<Mr) 

A  V  (crav) 

A  L  («rAU) 

A  V  (<7M-) 

A  L  (erAl?) 

A  V  (<ra„) 
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3  -2 

2.8(0.6) 

2.6(1.3) 

4.1(1.8) 

9.1(3.1) 

6.3(0.7) 

9.6(1.4) 

2  -  3 

2.6(0.4) 

3.2(0.9) 

6.3(1.6) 

3.7(1.0) 

5.9(0.4) 

6.7(0.6) 

From  Table  14,  after  using  the  FISHER-BEHERENS  test  to  compare  AU 
from  MINI  RANGER-LORAN  (TD)  against  AU  from  MINI 
RANGER-LORAN  (DISPLAY)  and  AV  against  AV  we  see  that  these  are  sta- 
tistically different  at  a  95%  confidence  interval  in  2  out  of  4  cases. 

Table  15.  SHIP  VELOCITY  COMPONENTS  FOR  MINI  RANGER,  LORAN 
(TD),  LORAN  (DISPLAY),  BOTTOM  TRACKING  FROM  POINT  3  TO 
POINT  4  AT  3  MINUTE  TIME  INTERVAL 


TIME 

MINI  RANGER 

LORAN  (TD) 

LORAN  (DIS- 
PLAY) 

BOTTOM 
TRACKING 

U 

(cm  sec) 

V 

(cm  sec) 

L 

(cm  sec) 

V 

(cm  sec) 

L 
(cm  sec) 

V 

(cm  sec) 

L 

(cm  sec) 

V 

(cm  sec) 

17  22  17 

-9.4 

286.1 

-12.1 

283.2 

-10.2 

275." 

0.5 

^o  ->  -> 

17  25  17 

-11.7 

288.0 

-12.8 

284.2 

-10.6 

288.9 

-1.4 

285.9 

17  28  17 

-12.5 

289.3 

-9.8 

285.3 

-10.9 

28". 5 

-4.1 

28".0 

17  31  17 

-4.9 

295.7 

-8.7 

295.9 

-15.1 

293.5 

1.5 

290.6 

17  34  17 

-8.7 

293.5 

-11.7 

292.8 

-".1 

293.0 

-1.9 

288.4 

17  37  19 

-19.4 

293.5 

-18.3 

286.7 

-23.6 

274.5 

-11.8 

290.1 

1 7  40  1 8 

-21.6 

29". 3 

-20.0 

296.2 

-20.0 

307.0 

-15.1 

293.2 

17  43  1 8 

-16.2 

288.6 

-15.9 

288.8 

-7.9 

288.8 

-9.2 

2S6.6 

1"  46  18 

-24.9 

292.2 

-26.4 

290.5 

-29.5 

299.3 

-20.2 

28".9 

17  49  18 

-37.5 

290.6 

-42.0 

284.7 

-41.6 

278.8 

-33.7 

286.8 

17  52  IS 

-32.2 

2S8.7 

-34.5 

291.0 

-41.3 

28". 5 

-2"  " 

286.4 

17  55  17 

-2S.0 

286.0 

-29.9 

281.1 

-23.4 

2S_\0 

-21.4 

2S2.5 

17  58  18 

-22  5 

281.3 

-13.6 

286.1 

-16.6 

24". 7 

-15.3 

278.6 

18  01  18 

-22.3 

281.0 

-23.4 

2"6.0 

-25.7 

280.0 

-18.5 

276.7 

18  04  18 

-25.3 

275.4 

-26.1 

275.1 

-20.0 

283.2 

-17.7 

273.9 

18  0"  IS 

-23.4 

274.6 

-24.9 

->--»-> 

-28. 3 

266.3 

-18.2 

271.6 
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18  10  18 

-22.3 

271.0 

-19.2 

273.1 

-22.3 

258.7 

-14.8 

268.9 

18  13  18 

-28.7 

267.2 

-30.2 

267.5 

-23.8 

285.7 

-23.1 

266.0 

18  16  18 

-29.8 

262.8 

-37.4 

257.3 

-37.7 

262.5 

-21.4 

265.1 

18  19  IS 

-30.9 

264.3 

-20.7 

269.9 

-20.0 

268.3 

-17.1 

262.9 

18  22  47 

-25.8 

261.2 

-33.0 

257.9 

-32.3 

256.5 

- 

- 

1 8  25  47 

-31.3 

256.1 

-37.8 

253.8 

-41.2 

247.6 

- 

- 

18  28  47 

-26.8 

266.0 

-29.1 

262.2 

-26.1 

262.7 

- 

- 

18  31  47 

-29.1 

263.0 

-25.7 

264.4 

-28.0 

251.7 

- 

- 

18  34  47 

-29.1 

264.7 

-30.2 

263.2 

-24.5 

281.2 

- 

- 

18  37  47 

-27.5 

263.9 

-27.2 

262.0 

-24.9 

245.9 

- 

- 

18  40  47 

-33.6 

261.1 

-35.8 

257.6 

-32.8 

264.4 

- 

- 

18  43  46 

-39.4 

256.5 

-38.2 

261.5 

-32.5 

251.1 

- 

- 

18  46  46 

-32.9 

259.8 

-34.4 

256.8 

-45.4 

269.8 

- 

- 

18  49  46 

-41.9 

257.5 

-40.8 

258.9 

-41.9 

258.0 

- 

- 

18  52  46 

-43.3 

256.0 

-43.3 

257.9 

-44.9 

262.8 

- 

- 

18  55  4" 

-37.7 

257.9 

-36.2 

255.3 

-33.2 

248.9 

- 

- 

18  58  4" 

-32.5 

255.4 

-31.0 

259.6 

-32.5 

262.2 

- 

- 

19  01  47 

-33.9 

252.4 

-32.8 

251.9 

-20.7 

231.4 

- 

- 

19  04  47 

-32.4 

255.7 

-32.4 

253.1 

-49.4 

274.3 

- 

- 

19  07  46 

-32.1 

253.1 

-32.1 

253.1 

-28.7 

237.3 

- 

- 

19  10  46 

-43.5 

253.2 

-41.2 

255.3 

-41.6 

273.4 

- 

- 

19  13  47 

-47.8 

256.5 

-48.1 

253.2 

-49.3 

252.9 

- 

- 

19  16  46 

-53.7 

255.4 

-56.4 

251.4 

-41.2 

232.0 

- 

- 

19  19  46 

-57.1 

256.5 

-56.3 

253.0 

-65.8 

273.3 

- 

- 

19  22  47 

-46.2 

259.8 

-51.9 

254.7 

-45.5 

242.7 

- 

19  25  46 

-42.0 

263.6 

-45.8 

259.4 

-53.3 

257.7 

- 

- 

19  2S  46 

-4". 7 

268.0 

-49.6 

264.9 

-49.6 

274.2 

- 

- 

19  31  46 

-50.9 

2"0.9 

-57.6 

263.3 

-49.7 

262.8 

- 

- 

19  34  47 

-59.5 

267.2 

-60.6 

262.5 

-65.9 

258.0 

- 

- 

19  37  46 

-54.9 

268.8 

-59.0 

262.9 

-53.4 

252.9 

- 

- 
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19  40  47 

-52.7 

272.5 

-51.9 

268.0 

-65.9 

283.8 

- 

19  43  47 

-60.4 

276.0 

-66.0 

270.3 

-66.0 

273.6 

- 

- 

19  46  46 

-57.9 

279.1 

-64.4 

272.7 

-61.7 

273.4 

- 

- 

19  49  46 

-45.2 

273.9 

-52.7 

268.4 

-52.3 

264.4 

- 

- 

19  52  47 

-39.1 

276.8 

-36.9 

266.9 

-34.2 

266.9 

- 

- 

19  55  47 

-21.1 

278.8 

-23.0 

275.8 

-28.6 

273.6 

- 

- 

19  58  46 

-21.6 

281.4 

-25.0 

276.6 

-16.3 

273.3 

- 

- 

20  01  47 

-21.8 

276.8 

-26.7 

274.2 

-33.1 

278.7 

- 

- 

20  04  47 

-19.5 

277.7 

-22.2 

271.3 

-20.7 

268.3 

- 

- 

20  07  47 

-26.8 

279.4 

-30.5 

269.9 

-24.5 

268.5 

- 

- 

20  10  46 

-23.4 

274.2 

-24.6 

274.9 

-25.0 

268.5 

- 

- 

Table   16.  COMPARISON  BETWEEN  MINI  RANGER  AND  LORAN  (TD)  AND 
THE  MINI  RANGER  AND  LORAN  (DISPLAY)  FOR  LEG  3  -  4. 


LEG 

MINI  RANGER  -  LORAN 
(TD) 

MINI  RANGER  -LORAN 
(DISPLAY) 

MINI  RANGER  -  BOT- 
TOM TRACKING 

A  U  (ff4U) 

ivy 

A  V  (<7S<) 

A  V  {cty) 

A  U  (aAa) 

A  V  (cM.) 

3  -4 

2.8(0.3) 

3.5(0.3) 

5.1(0.5) 

8.6(0.8) 

7.1(0.5) 

3.1(0.3) 

Table  15  again  yields  similar  results  for  leg  3  -  4.  From  Table  16  it  is  again 
obvious  that  the  results  which  comes  from  the  LORAN  (TD)  and  the  LORAN 
(DISPLAY)  are  different. 

In  Table  17   are  the  average  velocity  components  for  the  different  positioning 

systems   and  for  the  bottom  tracking  for  each  leg.  From  this  table  we  can  see: 

•  Data  averaged  over  times  of  about  18  minutes  (minimum)  show  no  statis- 
tically significant  difference  between  MINI  RANGER,  LORAN  (TD)  and 
LORAN  (DISPLAY).  Any  one  of  these  positioning  systems  can  be  used  for 
these  averaged  times.  This  is  due  to  the  fact  that  the  variations  of  the  ship 
course  and  speed  in  the  averaged  period  of  time  completely  swamp  the  dif- 
ferent precision  which  the  different  navigation  systems  have  in  their  calcu- 
lation of  the  U  and  V  components  of  the  ship.  The  minimum  time  interval 
from  which  the  averaged  data  over  the  different  navigation  systems  become 
the  same  statistically  depends  upon  the  ship's  variations  in  course  and  speed. 
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•  The  results  from  Bottom  Tracking  are  clearly  systematically  different  from 
those  of  the  others  systems.  They  do,  however,  show  the  same  level  of 
internal  consistency  as  both  the  MINI  RANGER  and  the  LORAN  (TD) 
results. 

•  The  much  larger  standard  deviations  on  the  LORAN  (DISPLAY)  results 
are  indicative  of  the  much  greater  variability  in  computed  U  and  V  compo- 
nents when  using  this  system. 


Table  17.    AVERAGE    VELOCITY    COMPONENTS    FOR       MINI    RANGER, 
LORAN  (TD),  LORAN  (DISPLAY),  BOTTOM  TRACKING 


MINT  RANGER 

LORAN  (TD) 

LORAN  (DIS- 
PLAY) 

BOTTOM  TRACK- 
ING 

<-  (frv) 

V  (fry) 

U  (*v) 

V  (C7V) 

u  K.) 

V  (ffy) 

V  (<rv) 

V(<Ty) 

1  - 

2 

-117.7(3.4) 

-279.9(1.2) 

-117.2(3.3) 

-276.7(2.1) 

-116.8(3.1) 

-275.6(3.2) 

■ 

- 

•> 
1 

146.1(3.7) 

256.3(1.9) 

144.7(3.8) 

252.7(1.4) 

143.3(3.2) 

250.6(4.1) 

- 

- 

1  - 
3 

-270.1(1.1) 

-122.3(1.1) 

-268.9(1.4) 

-121.2(1.3) 

-268.5(2.9) 

-121.4(3.2) 

- 

- 

3  - 
2 

276.8(1.1) 

-99.3(2.1) 

276.2(2.2) 

-97.5(1.5) 

274.5(3.1) 

-96.6(4.9) 

270.5(1.6) 

-108.9(1.8; 

-> 
3 

-264.7(1.4) 

119.1(1.1) 

-263.3(1.4) 

117.7(1.7) 

-264.1(3.0) 

118.7(2.1) 

-258.8(1.6) 

125.8(0.8) 

3  - 

4 

-32.5(1.8) 

271.1(1.7) 

-33.8(1.9) 

269.1(1.7) 

-33.6(2.1) 

268.7(2.1) 

- 

- 

3  - 
4A 

-21.6(1.9) 

283.3(2.4) 

-21.8(2.1) 

281.8(2.3) 

-21.8(2.3) 

282.1(2.7) 

-14.5(2.1) 

280.6(2.1) 

In  Appendix  E  there  are  12  figures  which  show  the  average  current  profiles 
from  point  1  -  3,  3  -  2,  2  -  3,  3  -  4  with  the  three  different  navigation  data  from 
MINI  RANGER,  LORAN  (TD).  LORAN  (DISPLAY).  These  plots  show  us 
the  same  thing  as  the  Table  17  i.e.  no  significant  difference  between  the  naviga- 
tion systems  when  we  average  over  times  of  about  20  -  30  minutes. 

In  Appendix  F  there  are  6  plots  for  points  1  -  2.  2  -  1  but  not  averaged  for 
the  whole  length  (just  plots  for  specific  times  between  the  points).  These  plots 
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show,  as  expected,  that  the  U  and  V  components  of  the  current  are  different  for 
different  navigation  systems. 
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• 


• 


VII.     CONCLUSIONS  AND  RECOMMENDATIONS 

From  this  thesis  we  can  conclude  that: 

Because  of  the  variability  in  the  LORAN  (DISPLAY)  results  and  the  need 
for  a  reasonable  measure  of  statistical  strength,  it  is  not  recommended  that 
LORAN  (DISPLAY)  be  used  for  giving  averaged  results  of  under  20  min- 
utes. Indeed,  it  is  preferable  to  use  LORAN  (TD)  rather  than  LORAN 
(DISPLAY)  whenever  possible.  Furthemore,  LORAN  (DISPLAY)  results 
not  only  show  a  much  greater  variability  than  the  LORAN  (TD)  results,  but 
also  at  a  95%  confidence  interval  consistently  fail  the  statistical  hypothesis 
that  they  come  from  the  same  population. 

MINI  RANGER  results  as  expected  are  very  consistent.  It  is  felt  that  good 
ADCP  current  estimates  could  be  derived  using  MINI  RANGER  observa- 
tions for  periods  as  short  as  3  minutes. 

When  averaging  for  20  -  30  minutes  all  the  navigation  systems  give  the  same 
results. 


•  Without  averaging  we  recommended  that  the  MINI  RANGER  be  used 
whenever  possible.  If  this  system  is  not  available,  LORAN  (TD)  should  be 
the  next  choice. 

•  The  theoritical  investigation  revealed  that  when  time  averaging  over  20 
minutes  or  longer,  the  choice  of  the  filter  to  be  used  for  the  LORAN  system 
was  not  critical.  The  cruise  results  appear  to  support  this  contention. 

•  Bottom  Tracking  results  are  systematically  biased,  but  have  the  same 
standard  deviation  as  the  MINI  RANGER  results. 

In  looking  to  the  future  work  on  this  subject  it  is  recommended  that: 


• 


Further  work  is  required  to  understand  and  resolve  the  Bottom  Tracking 
bias  problem. 


•  If  this  cruise  is  done  again  we  suggest  checking  the  LORAN  filters  by  mak- 
ing one  long  run,  dividing  it  into  30  minute  segments,  and  using  a  different 
filter  on  each  segment. 
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APPENDIX  A.     STATISTICAL  ANALYSIS  FOR  THE  LC  408  DATA 


The  provided  numbers  are  in   nanoseconds  . 

(1)     TIME  CONSTANT  0  SECONDS  (NO  FILTERING) 


STATION 

N 

MEAN 

WHISKY 

1440 

16312005. 09 

X-RAY 

1440 

27511618.59 

YANKEE 

1440 

42735925. 76 

(2)     TIME  CONST 

STATION 

N 

MEAN 

WHISKY 

1440 

16312033. 94 

X-RAY 

1440 

27511616. 71 

YANKEE 

1440 

42735928.43 

(3)     TIME  CONST 

STATION 

N 

MEAN 

WHISKY 

1440 

16311981. 04 

X-RAY 

1440 

27511616.56 

YANKEE 

1440 

42735918.06 

STANDARD 
DEVIATION 

46. 85663572 
22.27564796 
24.94991412 


STANDARD 
DEVIATION 

51.  75565604 
19. 82896815 
23. 39278554 


MINIMUM 
VALUE 

16311760.00 
27511540.00 
42735810.00 


MINIMUM 
VALUE 

16311860. 00 
27511560. 00 
42735860. 00 


MAXIMUM 
VALUE 

16312180.00 
27511720. 00 
42736030. 00 


MAXIMUM 
VALUE 

16312180.00 
27511680. 00 
42736020. 00 


STD  ERROR 
OF  MEAN 

1.23478077 
0.58701487 
0.65748797 


STD  ERROR 
OF  MEAN 

1.  36388129 
0. 52253919 
0. 61645403 


TIME  CONSTANT  10  SECONDS  (SLOWER  RESPONSE) 


STANDARD 
DEVIATION 

49. 40674807 
13. 76653231 
18. 02041879 


(4)     TIME 
APPLICATIONS; 


CONSTANT 


W 


STATION 


WHISKY 

X-RAY 

YANKEE 


N 


MEAN 


1417  16311990.54 
1417  27511618.80 
1417   42735915.31 


STANDARD 
DEVIATION 

46. 87615596 
12. 99066592 
16. 74596778 


MINIMUM 
VALUE 

16311790. 00 
27511580. 00 
42735850. 00 


MAXIMUM 
VALUE 

16312110. 00 
27511680. 00 
42735990. 00 


STD  ERROR 
OF  MEAN 

1. 30198213 
0. 36277998 
0.47487973 


SECONDS        (MONITOR 


MINIMUM 
VALUE 

16311800. 00 
27511570. 00 
42735850. 00 


MAXIMUM 
VALUE 

16312130. 00 
27511680. 00 
42735980.00 


STD  ERROR 
OF  MEAN 

1.  24528015 
0. 34510122 
0.44486202 
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APPENDIX  B.     TABLES  FOR  STANDARD  DEVIATIONS  AND 
COVARIANCES  OF  VELOCITY  COMPONENTS  FOR  LORAN  C  (LC  408) 

Table  18.       TIME  INTERVAL  BETWEEN  POSITIONS  3  MINUTES 


FILTER 

cv 

<*U 

cv,u 

0  sec 

1 1.6  cm  sec 

6.8  cm/sec 

-61.184  cm2/ sec2 

5  sec 

10.8  cm  sec 

6.1  cm /sec 

-51.047  cm2 1  sec2 

10  sec 

8.2  cm/sec 

4.3  cm/sec 

-27.447  cm2/  sec2 

20  sec 

7.7  cm-sec 

4.1  cm/sec 

-24.029  cm2/  sec2 

Table  19.       TIME  INTERVAL  BETWEEN  POSITIONS  6  MINUTES 


FILTER 

GV 

oL 

°v,v 

0  sec 

5.8  cm /sec 

3.4  cm -see 

-15.296  cm2/  sec2 

5  sec 

5.4  cm/sec 

3.0  cm/sec 

-12.761  cm2/  sec2 

10  sec 

4.1  cm  sec 

2.1  cm /sec 

-6.861  cm2/  sec2 
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20  sec 

3.8  cm/sec 

2.0  cm/sec 

-6.007  cm2/  sec2 

Table  20.       TIME  INTERVAL  BETWEEN  POSITIONS  12  MINUTES 


FILTER 

Oy 

<>u 

°V,U 

0  sec 

2.9  cm,  sec 

1.7  cm  sec 

-3.824c/??2/  sec2 

5  sec 

2.7  cm  sec 

1 .5  cm/sec 

—3.190c/??2/  sec2 

10  sec 

2.1  cm  sec 

1.1  cm  sec 

-1.715c/??2/  sec2 

20  sec 

1.9  cm  sec 

1.0  cm,sec 

—  1.502c/??2/  sec2 

Table  21.       TIME  INTERVAL  BETWEEN  POSITIONS  20  MINUTES 


FILTER 

GV 

cu 

0v,u 

0  sec 

1 .7  cm  'sec 

1.0  cm/ sec 

-1.376c//?2/  sec2 

5  sec 

1.6  cm,: sec 

0.9  cm  sec 

-1.148c//?2/  sec2 
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10  sec 

1 .2  cm/sec 

0.7  cm/sec 

— 0.61 1cm2/  sec2 

20  sec 

1.1  cm/sec 

0.6  cm/sec 

-0.540cm2/  sec2 

Table  22.       TIME  INTERVAL  BETWEEN  POSITIONS  40  MINUTES 


FILTER 

GV 

°v 

0v,u 

0  sec 

0.9  cm  sec 

0.5  cm  sec 

-0.344cm2/  sec2 

5  sec 

0.8  cm  sec 

0.4  cm  sec 

-0.287cm2/ sec2 

10  sec 

0.6  cm  sec 

0.3  cm  sec 

-  0.154cm2/  sec2 

20  sec 

0.5  cm/sec 

0.3  cm  sec 

-0.135cm2/ sec2 

50 


APPENDIX  C.     PROGRAM  DRIVLR  FORTRAN 

A.  GENERAL  REMARKS 

This  program  calculate  the  hyperbolic  coordinates  of  a  ship  expressed  in 
LORAN  time  differences  from  geographic  potisions  forward  computation  ,  and 
the  inverse  computation. 

Forward  computation: 

•  Input  parameters  -  ALAT,  ALON,  L,  DEL,  ISYS,  IB,  IVEL. 

•  Output  parameters  -  TD,  BETA,  BLEN,  BLEM,  IER. 

Inverse  computation: 

•  Input  parameters  -  ALAT.  ALON,  L,  DEL,  TD.  ISYS,  ISENT.  ISW.  IB. 
IVEL. 

•  Output  parameters  -  ALAT,  ALON,  BETA,  BLEN,  BLEM,  ISENT.  IER. 
IT. 

B.  DESCRIPTION  OF  PARAMETERS: 

ALAT  -  Latitude  array  in  seconds  of  arc: 
ALAT(l)  -  Latitude  of  master  station. 
ALAT(2)  -  Latitude  of  slave  1  station. 
ALAT(3)  -  Latitude  of  master  station. 
ALAT(4)  -  Latitude  of  slave  2  station. 
ALAT(5)  -  Latitude  of  ship  position. 
ALON(l)  -  Longitude  of:  master  station. 
ALON(2)  -  Longitude  of  slave  1  station. 
ALON(3)  -  Longitude  of  master  station. 
ALON(4)  -  Longitude  of  slave  2  station. 
ALON(5)  -  Longitude  of  ship  position. 

L  -  Spheroid  code  array. 
L(l)  -  Spheroid  code  array  for  first  rate: 
L(l)  =  1      WGS  1984. 
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L(l) 

=  2 

Bessel  spheroid. 

L(l) 

=  3 

Clarke  1858  spheroid. 

L(l) 

=  4 

Clarke  1866  spheroid. 

L(l) 

=  5 

Clarke  1880  spheroid. 

L(l) 

=  6 

Everest  spheroid. 

L(l) 

=  7 

Fischer  spheroid. 

L(l) 

=  8 

International  spheroid. 

L(l) 

=  9 

WGS  1972   spheroid. 

-(2)- 

Spheroid  code  for  second  rate: 

L(2) 

=  1 

WGS  1984. 

L(2) 

Bessel  spheroid. 

L(2) 

=  3 

Clarke  1858  spheroid. 

L(2) 

=  4 

Clarke  1866  spheroid. 

L(2) 

=  5 

Clarke  1880  spheroid. 

L(2) 

=  6 

Everest   spheroid. 

L(2) 

=  7 

Fischer  spheroid. 

L(2) 

=  8 

International  spheroid. 

L(2) 

=  9 

WGS  1972   spheroid. 

DEL  -  Coding  delay  array,  in  microseconds: 
DEL(l)        Coding  delay  for  first  rate. 
DEL(2)        Coding  delay  for  second  rate. 

TD  -  Time  difference  array,  in  microseconds: 
TD(1)         Time  difference  for  first  rate. 
TD(2)         Time  difference  for  second  rate. 

I  SYS  -  Loran  system  array: 

ISYS(l)       Loran  svstem  for  first  rate. 
ISYS(1)'=A:  Loran  A 
ISYS(l)  =  c:  Loran  C 

ISYS(2)       Loran  svstem  for  second  rate. 
ISYS(2)'=A:  Loran  A 
ISYS(2)  =  c:  Loran  C 
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BETA  -  Baseline  delay  array: 
BETA(l)       Baseline  delay  for  first  rate. 
BETA(2)       Baseline  array  for  second  rate. 

BLEN  -  Minimum  time  difference  arrav,  in  microseconds: 
BLEN(l)       Minimum  time  difference  for  first  rate. 
BLEN(2)       Minimum  time  difference  for  second  rate. 

BLEM  -  Maximum  time  difference  array,  in  microseconds: 
BLEM(l)       Maximum  time  difference  for  first  rate. 
BLEM(2)       Maximum  time  difference  for  second  rate. 

ISENT  -  Sentinel  for  DR  position: 
ISENT  =  0       DR  position  not  given. 
ISENT  =  1        DR  position  given. 

ISW        Operation  code: 
ISW  =  1  Forward  computation. 

ISW  =  2         Inverse  computation. 

IB  -  Switch  for  baseline  computation: 
IB  =  0  Baseline  computation  off. 

IB  =  1  Baseline  computation  on. 

IER  -  Error  code: 

I ER  =  1  Invalid  velocity  of  propagation  code. 

IER  =  2  Invalid  spheroid  code. 

IER  =  3  Invalid  configuration. 

IER  =  4  Invalid  operation  code. 

IER  =  5  Invalid  readings. 

IER  =  6  Invalid  DR  position. 

IT  -  Iteration  code. 

IVEL  -  Velocity  of  propagation  code  array: 

IVEL(l)  -  Velocity  of  propagation  code  for  first  rate: 
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•  0  =  299.6929      kilometers  per  microsecond 

•  1  =  299.691 1624  kilometers  per  microsecond 

•  2  =  299.708       kilometers  per  microsecond 

IVEL(2)  -  VElocity  of  propagation  code  for  second  rate: 

•  0  =  299.6929      kilometers  per  microsecond 

•  1  =  299.691 1624  kilometers  per  microsecond 

•  2  =  299.708       kilometers  per  microsecond 

C  *  PROGRAM  DRIVLR  FORTRAN  SOURCE  CODE  * 

C 

C  DESCRIPTION  OF  PARAMETERS 

C         DATA  DECLARATION 

IMPLICIT  REAL*8(A-H,0-Z) 

REAL*8  ALAT(5),ALON(5),DEL(2),TD(2),BETA(2),BLEN(2),BLEM(2),LN 
C     CHARACTER--'"  1  ISYS(2) 

INTEGER  H,M,L(2),IVEL(2),ISYS(2) 

ALAT( 1)=3600.  0D0*39.  0D0  +  33. 0D0*60. 0D0  +  6. 740D0 

ALAT(2)=3600. 0D0-38.  0D0  +  46. 0D0*60. 0D0  +  57.  110D0 

ALAT(3)=3600. 0D0*39. 0D0  +  33. 0D0*60. 0D0  +  6. 740D0 

ALAT(4)=3600.0D0*35.0D0  +  19. 0D0*60. 0D0  +  18.  3050D0 

ALAT(5)=3600.0D0*36. 0D0  +  45. 0D0*60. 0D0  +  0.  0D0 

ALON(l)=3600.0D0*118. 0D0  +  49. 0D0*60. 0D0  +  55. 8160D0 

ALON(2)=3600. 0D0*122. 0D0  +  29. 0D0*60. 0D0  +  43. 9750D0 

ALON(3)=3600. 0D0*118. 0D0  +  49. ODO*60. 0D0  +  55. 8160D0 

ALON(4)=3600. ODO-114. 0D0  +  48. 0D0*60. 0D0  +  16. 8810D0 

ALON(5)=3600. 0D0*121. 0D0  +  55. 0D0*60. 0D0  +  0.  0D0 

L(l)  =  1 

L(2)  =  1 

DEL(l)  =  27000. 0D0 

DEL(2)  =  40000. 0D0 

ISYS(1)=1 

ISYS(2)=1 

ISENT=1 

ISW=2 

IVEL(1)=1 

IVEL(2)=1 

IB=1 

CALL  LORAN  ( ALAT, ALON ,L,DEL,TD, I SYS ,BETA,BLEN,BLEM, I SENT, ISW, IB , 
*IER,IT,IVEL) 

DO  200  1=1,914 

READ(7,10)  H,M,S,TD(1),TD(2) 
10   F0RMAT(I2,1X,I2,1X,F5. 2,21X,F8.  2,1X,F8.  2) 

IER=0 

CALL  LORAN  (ALAT, ALON, L, DEL, TD, ISYS , BETA, BLEN,BLEM, ISENT, ISW, IB , 
*IER,IT,IVEL) 

WRITE(9,20)  I,H,M,S,ALAT(5),ALON(5),IT,IER 
20  FORMAT(I3,1X,I2,1X,I2,1X,F5.2,2X,F10.2,2X,F10.2,2X,I2,2X,I1) 

F=ALAT(5)/3600.  0D0 

LN=-ALON(5)/3600.  0D0 
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WRITE(10,30)  I,H,M,S,F,LN 
30  F0RMAT(I3,1X,I2,1X,I2,1X,F5.  2,2X,F10. 7 ,2X,F12.  7) 

READ(8,40)  F1,LNI 
40  FORMAT(9X,F10.  7,1X,F12.  7) 

DF=F-F1 

DL=LN-LN1 

WRITE(6,*)  DF,DL,LN1 
200  CONTINUE 

STOP 

END 

SUBROUTINE  LORAN  (  ALAT, ALON,L,DEL,TD, ISYS,BETA,BLEN,BLEM, 
*ISENT,ISW,IB,IER,IT,IVEL) 
C 

C        DATA  DECLARATION 
C 

IMPLICIT  REAL*8(A-H,0-Z) 

INTEGER  SOUTH, EAST, WEST 
C 

C        STORAGE  ALLOCATION 
C 

DIMENSION  DEL(2),TD(2),BETA(2),BLEN(2),BLEM(2),ISYS(2),H(2),R(2), 
*IVEL(2),VELL(2),L(2) 

DIMENSION  VEL(3) 

DIMENSION  TEMP1(4),TEMP2(4) 

DIMENSION  ALAT(5),ALON(5) 

DIMENSION  A(9),B(9) 

DATA  A/  6378137.  0D0  ,  6377397.  1550D0  ,  6378293. 6450D0  , 
*  6378206. 4000D0  ,  6378249. 1450D0, 6377276. 3450D0 ,6378166. 0000D0, 
*6378388.  0000D0 , 6378135. 0000D0/ 

DATA  B/6356752. 31420D0 , 6356078. 9628D0 , 6356617. 9376D0, 
*6356583.  8000D0 ,6356514.  8695D0 ,6356075. 4131D0 ,6356784. 2836D0, 
*6356911. 9461D0, 6356750. 5200D0/ 

DATA  Al/24.  0305D0/ , A2/-0. 40758D0/ ,A3/0. 00346776D0/ 

DATA  Bl/0.  510483D0/,B2/-0. 011402D0/ ,B3/0. 001760D0/ 

DATA  ARC1/0.484813681110D-5/ 

DATA  LORANC/1/ 

DATA  FAZ/0. 0D0/ , IK/ 0/ , TEMP 1/4*0.  0D0/ ,TEMP2/4*0.  0D0/ 

DATA  ZERO/0. 0D0/ 

DATA  N0RTH,S0UTH,EAST,WEST/1HN,1HS,1HE,1HW/ 

DATA  MINUS/ 1H-/ 

DATA  VEL/299. 6929D0,299. 6911624D0 ,299.  708D0/ 

DATA  BAZ/0. ODO/,JER/0/ 

DATA  ALATS/0. 0D0/ ,ALONS/0. ODO/,N/0/ 
C 

C        INITIALIZE  ERROR  CODE 
C 

IER=0 
C 

C        TEST  SWITCH  FOR  BASELINE  COMPUTATION 
C 

IF  (IB)  5,150,5 
C 

C        BASELINE  COMPUTATION 
C 
C 
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C        TEST  FOR  VALID  VFLOCITY  OF  PROPAGATION  CODE  FOR  BOTH  CHAINS 
C 

5   DO  10  1=1,2 

IF  (IVEL(I).GT.  2)  GO  TO  15 
10  CONTINUE 
GO  TO  20 
15   IER=1 

GO  TO  320 
C 

C        DETERMINE  VELOCITY  OF  PROPAGATION  FOR  BOTH  CHAINS 
C 

20  DO  25  1=1,2 

IVEL(I)=IVEL(I)+1 
KK=IVEL(I) 
25   VELL(I)=VEL(KK) 
A4A=VELL(1) 
A4B=VELL(2) 
C 

C        TEST  FOR  CHANGE  IN  SPHEROID  CODE,  FOR  VALID  SPHEROID  CJODE 
C        AND  DETERMINE  APPROPRIATE  SPHEROID  COSTANTS 
C 

DO  30  J=l,2 

IF  (L(J).GE.  LAND.  L(J).LE.9)  GO  TO  30 
GO  TO  55 
30  CONTINUE 

IF  (L(l).EQ.  L(2))  GO  TO  35 
GO  TO  55 
35   K=L(1) 

IF  (IK-K)  40,60,40 
40  DO  50  1=1,9 

IF  (I-K)  50,45,50 
45   A5=A(I) 
A6=B(I) 
A55=A5*A5 

A10=(A55-A6*A6)/A55 
IK=K 

GO  TO  60 
50  CONTINUE 
55   IER=2 

GO  TO  320 
C 

C        TEST  FOR  CHANGE  IN  FIXED  STATIONS 
C 

60   DO  80  1=1,4 

IF  (TEMPl(I)-ALAT(I))  70,65,70 
65  IF  (TEMP2(I)-ALON(I))  70,80,70 
70   DO  75  J=l,4 

TEMP1(J)=ALAT(J) 
75   TEMP2(J)=ALON(J) 

GO  TO  85 
80   CONTINUE 
GO  TO  145 
C 

C        TEST  FOR  TRIAD  OR  TETRAD  CONFIGURATION  AND  ACTIVATE 
C        APPROPRIATE  SENTINEL 
C 
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85  A12=ALAT(1)-ALAT(3)+AL0N(1)-AL0N(3) 
A12=DABS(A12) 

A13=ALAT(2)-ALAT(4)+AL0N(2)-AL0N(4) 
A13=DABS(A13) 

IF  (A12-0. 001D0)  90,90,105 
90  IF  (A13-0. 001D0)  95,95,100 
95   IER=3 

GO  TO  320 
100  All=-1. ODO 

GO  TO  110 
105  All=l.  ODO 
■£ 

C        TEST  FOR  VALID  OPERATION  CODE 
C 

110  IF  (ISW-1)  120,125,115 
115  IF  (ISW-2)  125,125,120 
120  IER=4 

GO  TO  320 
C 

C        COMPUTE  GEODETIC  DISTANCE  AND  SALT-WATER  RETARDATION  FROM 
C        MASTER  1  STATION  TO  SLAVE  1  STATION 
C 

125  CALL  SODIN  ( ALAT( 1) , ALON( 1) , ALAT( 2) , ALON( 2) ,L, A35 ,FAZ,BAZ , JER) 
AA10=A35/A4A 
M=l 
1=1 

GO  TO  295 
130  AA11=A47 
C 

C        COMPUTE  GEODETIC  DISTANCE  AND  SALT-WATER  RETARDATION  FROM 
C        MASTER  2  STATION  TO  SLAVE  2  STATION 
C 

CALL  SODIN  (ALAT(3),AL0N(3),ALAT(4),AL0N(4),L,A35,FAZ,BAZ,JER) 
CC10=A35/A4B 
M=2 
1=2 

GO  TO  295 
135  CC11=A47 
C 

C        COMPUTATION  OF  BASELINE  DATA 
BETA(1)=AA10+AA11 
BETA(2)=CC10+CC11 
DO  140  1=1,2 
BLEN(I)=DEL(I) 
140  BLEM(I)=2.  ODO*BETA(I)+DEL(I) 
C 

C        TURN  SWITCH  FOR  BASELINE  COMPUTATION  OFF 
C 

145  IB=0 
RETURN 
C 

C        TEST  SWITCH  FOR  FORWARD  OR  INVERSE  COMPUTATION 
C 

150  GO  TO  (155,190),ISW 
C 
C        FORWARD  COMPUTATION 
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c 
c 

C        COMPUTE  GEODETIC  DISTANCE  AND  SALT-WATER  RETARDATION  FROM 

C        SHIP  TO  MASTER  1  STATION 

C 

155  CALL  SODIN  ( ALAT(5) ,ALON(5) ,ALAT( 1) ,ALON( 1) ,L,A35 ,FAZ,BAZ, JER) 

A76=A35/A4A 

A72=A76 

M=l 

1=7 

GO  TO  295 
160  A77=A47 

A73=A77 
C 

C  COMPUTE  GEODETIC  DISTANCE  AND  SALT-WATER  RETARDATION  FROM 

C  SHIP  TO  SLAVE  1  STATION 

C 

CALL  SODIN  (ALAT(5),ALON(5),ALAT(2),ALON(2),L,A35,FAZ,BAZ,JER) 

A74=A35/A4A 

M=l 

1=8 

GO  TO  295 
165  A75=A47 
C 

C        COMPUTE  GEODETIC  DISTANCE  AND  SALT-WATER  RETARDATION  FROM 
C        SHIP  TO  SLAVE  2  STATION 
C 

CALL  SODIN  (ALAT(5),ALON(5),ALAT(4),ALON(4),L,A35,FAZ,BAZ,JER) 

A70=A35/A4B 

M=2 

1=9 

GO  TO  295 
170  A71=A47 
C 

C        COMPUTE  GEODETIC  AND  SALT-WATER  RETARDATION  FROM 
C        SHIP  TO  MASTER  2  STATION 
C 

IF  (All)  185,175,175 
175  CALL  SODIN  ( ALAT(5) ,ALON(5) ,ALAT(3) ,ALON( 3) ,L,A35 ,FAZ,BAZ, JER) 

A72=A35/A4B 

M=2 

1  =  10 

GO  TO  295 
180  A73=A47 
C 

C        COMPUTE  LORAN  TIME  DIFFERENCES  OF  FIX 
C 

185  TD( 1)=AA10+AA11+A74+A75-A76-A77+DEL( 1) 

TD(2)=CC10+CC11+A70+A71-A72-A73+DEL(2) 

RETURN 
C 

C        INVERSE  COMPUTATION 
C 
C 

C        STORE  DR  POSITION  OF  SHIP 
C 
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190  ALATDR=ALAT(5) 

ALONDR=ALON(5) 
C 

C        TEST  FOR  VALID  MICROSECOND  READINGS 
C 

DO  195  1=1,2 

IF  (TD(I).LE.  BLEM(I).AND.TD(I).GE.  BLEN(I))  GO  TO  195 

IER=5 

GO  TO  320 
195  CONTINUE 
C 

C        CONVERT  SHIPS  POSITION  TO  RADIANS 
C 

IF  (ISENT)  205,205,200 
200  ALATR=ALAT(5)*ARC1 

ALONR=ALON(5)*ARCl 
C 

C        INITIALIZE  IETRATION  COUNTER 
C 

205  IT=0 

GO  TO  215 
C 

C        TEST  ITERATION  COUNTER 
C 

210  IF  (IT.  LT.  30)  GO  TO  215 

IER=6 

GO  TO  320 
C 

C        COMPUTE  SINES  AND  COSINES  OF  SHIPS  POSITION 
C 

215  SINS=DSIN(ALATR) 

COSS=DCOS(ALATR) 
C 

C        COMPUTE  GEODETIC  DISTANCE,  SIN  AND  COSINE  OF  AZIMUTH,  AND 
C        SALT-WATER  RETARDATION  FROM  SHIP  TO  SLAVE  2  STATION 
C 

CALL  SODIN  (ALAT(5),ALON(5),ALAT(4),ALON(4),L,A35,FAZ,BAZ,JER) 

M=2 

1=3 

GO  TO  280 
220  C1=A35 

D1=C1/A4B 

C2=DSIN(FAZ) 

C3=DCOS(FAZ) 

C101=A47 
C 

C        COMPUTE  GEODETIC  DISTANCE,  SIN  AND  COSINE  OF  AZIMUTH,  AND 
C        SALT-WATER  RETARDATION  FROM  SHIP  TO  MASTER  2  STATION 
C 

CALL  SODIN  (ALAT(5),ALON(5),ALAT(3),ALON(3),L,A35,FAZ,BAZ,JER) 

M=2 

1=4 

GO  TO  280 
225  C4=A35 

D4=C4/A4B 

C5=DSIN(FAZ) 
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C6=DC0S(FAZ) 

C104=A47 
C 

C        COMPUTE  GEODETIC  DISTANCE,  SIN  AND  COSINE  OF  AZIMUTH,  AND 
C        SALT-WATER  RETARDATION  FROM  SHIP  TO  SLAVE  1  STATION 
C 

CALL  SODIN  (ALAT(5),ALON(5),ALAT(2),ALON(2),L,A35,FAZ,BAZ,JER) 

M=l 

1=5 

GO  TO  280 
230  C7=A35 

E7=C7/A4A 

C8=DSIN(FAZ) 

C9=DCOS(FAZ) 

C107=A47 

IF  (All)  245,235,235 
C 

C        COMPUTE  GEODETIC  DISTANCE,  SIN  AND  COSINE  OF  AZIMUTH,  AND 
C        SALT-WATER  RETARDATION  FROM  SHIP  TO  MASTER  1  STATION 
C        (FOR  TETRAD  CONFIGURATIONS  ONLY) 
C 

235  C10=C7 

C11=C8 

C12=C9 

CALL  SODIN  (ALAT(5),AL0N(5),ALAT(1),AL0N(1),L,A35,FAZ,BAZ,JER) 

M=l 

1=6 

GO  TO  280 
240  C7=A35 

D7=C7/A4A 

C8=DSIN(FAZ) 

C9=DCOS(FAZ) 

C108=A47 
C 

C        DETERMINATION  OF  TIME  DIFFERENCE  COSTANTS  AND  DIFFERENTIAL 
C        CORRECTIONS  IN  Y  AND  X 
C        (FOR  TETRAD  CINFIGURATIONS  ONLY) 
C 

C13=TD(2)-CC10-CC11-C101+C104-DEL(2) 

C17=C13*A4B 

C18=TD( 1) -AA10-AA11-C107+C108-DEL( 1) 

C22=C18*A4A 

C23=C1-C17 

C24=C4 

C25=C7+C22 

C26=C10 

C27=((C2-C5)*(C25-C26)*(C23-C24)*(C11-C8)) 

C29=((C2-C5)*(C12-C9)*(C3-C6)*(C8-C11)) 

C30=C27/C29 

C28=(C23-C24+C30-«(C3-C6))/(C5-C2) 

GO  TO  250 
C 

C        DETERMINATION  OF  TIME  DIFFERENCE  COSTANTS  AND  DIFFERENTIAL 
C        CORRECTIONS  IN  Y  AND  X 
C        (FOR  TRIAD  CINFIGURATIONS  ONLY) 
C 
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245  C13=TD(2)-CC10-CC11-C101+C104-DEL(2) 
C17=C13*A4B 

C18=TD(1)-AA10-AA11-C107+C104-DEL(1) 
C22=C18*A4A 

C23=C1-C17 

C24=C4 

C25=C7-C22 

C27=C2*(C25-C24)+C23*(C5-C8)+C8*C24-C5*C25 

C29=C2*(C6-C9)+C3*(C8-C5)+C5*C9-C8*C6 

C30=C27/C29 

C28=(C23-C24+C30*(C3-C6))/(C5-C2) 
C 

C        DETRMINATION  OF  DIFFERENTIAL  CORRECTIONS  IN  LATITUDE  AND 
C        LONGITUDE  AND  PRELIMINARY  GEOGRAPHIC  POSITION  OF  SHIP 
C 

250  C100=l. 0D0-A10*SINS*SINS 

C102=DSQRT(C100) 

C32=A5/C102 

C31=C32*((1.0D0-A10)/C100) 

C33=C30/C31 

C34=(-C28/(C32*COSS)) 

ALATR=ALATR+C33 

AL0NR=AL0NR+C34 

ALAT(5)=ALATR/ARC1 

ALON(5)=ALONR/ARCl 
C     WRITE(13,*)  ALAT(5),ALON(5),'ALAT  ALON  FM  LORAN' 
C 

C        TEST  OF  CONVERGENCE  OF  SOLUTION 
C 

IF  (All)  260,255,255 
255  H(1)=AA10+AA11+E7+C107-D7-C108+DEL(1) 

GO  TO  265 
260  H( 1)=AA10+AA11+E7+C107-D4-C104+DEL( 1) 
265  H(2)=CC10+CC11+D1+C101-D4-C104+DEL(2) 

DO  270  1=1,2 
270  R(I)=TD(I)-H(I) 

ANUM=DABS(R(1))+DABS(R(2)) 

IF  (ANUM. LT. 0. 0001D0)  GO  TO  275 

IT=IT+1 

GO  TO  210 
275  RETURN 
C 

C        CONVERSION  OF  AZIMUTHS  FROM  NORTH  TO  SOUTH  ORIENTATION 
C 

280  IF  (FAZ-648000. 0D0)  285,285,290 
285  FAZ=(FAZ+648000. 0D0)*ARC1 

GO  TO  295 
290  FAZ=(FAZ-648000. 0D0)*ARC1 
C 
C 

C        FORWARD  AND  INVERSE  COMPUTATION 
C        OF  SALT-WATER  RETARDATION  (FOR  LORAN  C  ONLY) 
C 

295  IF  (ISYS(M).EQ. LORANC)  GO  TO  300 

A47=0. 0D0 

GO  TO  315 
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300  A46=A35/1609.  349D0 

IF  (A46-100. 0D0)  305,310,310 
305  A47=(B1/A46)+B2+B3*A46 

GO  TO  315 
310  A47=(A1/A46)+A2+A3*A46 
C 

C        BRANCHING  FOR  QUANTITIES  TO  BE  USED  IN  SOLUTION 
C 

315  GO  TO  (130, 135, 220, 225, 230, 240, 160, 165, 170, 180), I 
C 

C        PRINT  ROUTINE  FOR  INVALID  DATA 
C 

320  IF  (IER.  LE.4)  GO  TO  340 

CALL  ANGLE  (N,ALONDR, ISGN,LOND,LONM,ALONS) 
CALL  ANGLE  (N,ALATDR, ISGM,LATD,LATM,ALATS) 
IF  (ISGM.EQ.  MINUS)  GO  TO  325 
KLA=NORTH 
GO  TO  330 
325  KLA=SOU.TH 

330  IF  (ISGN.EQ. MINUS)  GO  TO  335 
KLO=WEST 
GO  TO  340 
335  KLO=EAST 

340  GO  TO  (345, 350, 355, 360, 365, 370), IER 
345  PRINT  375 

RETURN 
350  PRINT  380 

RETURN 
355  PRINT  385 

RETURN 
360  PRINT  390 

RETURN 
365  PRINT  395,LATD,LATM,KLA,LOND,LONM,KLO,TD(l),TD(2) 

RETURN 
370  PRINT  400,LATD,LATM,KLA,LOND,LONM,KLO,TD(1),TD(2) 
RETURN 
C 

C     FORMAT  STATEMENTS 
C 

375  F0RMATUH1,  IX,'  INVALID  VELOCITY  OF  PROPAGATION  CODE,  DATA  NOT  PROC 

*ESSED' ) 
380  F0RMAT(1H1, IX, 'INVALID  SPHEROID  CODE, DATA  NOT  PROCESSED') 
385  F0RMAT(1H1,1X,' INVALID  CONFIGURATION,  DATA  NOT  PROCESSED') 
390  F0RMAT(1H1, IX,' INVALID  OPERATION  CODE,  DATA  NOT  PROCESSED') 
395  F0RMAT(18X,2I3,A1,6X,I4,I3,A1,16X,F9.  2,13X,F9.  2,3X, 'INVALID  READIN 

*GS '  ) 
400  F0RMAT(18X,2I3,A1,6X,I4,I3,A1,16X,F9.  2,13X,F9.  2 ,3X, ' INVALID  DR  POS 
*ITION') 
END 


SUBROUTINE  SODIN  ( ALAT1 ,AL0N1 ,ALAT2,ALON2,L,S,FAZ,BAZ, IER) 
IMPLICIT  REAL*8(A-H,0-Z) 
DIMENSION  X(2),E(2),F(2),L(2) 
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DATA  PI/3. 14159265359D0/,ARCl/0.484813681110D-5/ 

DATA  A,B/6378137. ODO, 6356752.  3142D0/ 
C 

C     INITIALIZE  ERROR  CODE  AND  DETERMINE  COSTANTS 
C 

A14=l. 0D0-(B/A) 

A56=A14*A14 

A50=l. 0D0+A14+A56 

A51=A50-1.  ODO 

A52=A56/2.  ODO 

A53=-A51/2.  ODO 

A54=A56/16.0D0 

A55=A56/8.0D0 

A57=A56*1.  25 ODO 

A58=A56/4.  ODO 
C 

C     TEST  FOR  IDENTICAL  STATIONS 
C 

30  IF(ALAT1-ALAT2)  45,35,45 
35  IF(AL0N1-AL0N2)  45,40,45 
40   IER=2 

GO  TO  255 
C 

C     COMPUTATION  OF  PARAMETRIC  LATITUDE  OF  STATIONS 
C 

45   X1=ALAT1*ARC1 

Y1=AL0N1*ARC1 

X2=ALAT2 » ARC  1 

Y2=AL0N2*ARC1 

X(1)=X1 

X(2)=X2 

DO  50  1=1,2 

BETA=DATAN((1.0D0-A14)*DSIN(X(I))/DCOS(X(I))) 

E(I)=DSIN(BETA) 
50   F(I)=DCOS(BETA) 
C 

C     COMPUTATION  OF  SPHERICAL  DIFFERENCE  OF  LONGITUDES 
C 

A59=-Y2 

A60=-Y1 

C35=A59-A60 

C36=DABS(C35) 

IF(C36-PI)  60,55,55 
55   A16=2.0D0*PI-C36 

GO  TO  65 
60  A16=C36 
65   IF(A16)  75,70,75 
70   A16=0.5D-7 
C 

C     COMPUTATION  OF  GEODETIC  DISTANCE 
C 

75   A17=DSIN(A16) 

A18=DC0S(A16) 

A19=E(1)»E(2) 

A20=F(1)*F(2) 

A21=A19+A20*A18 
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A80=A17*F(2) 

A81=A80*A80 

A82=E(2)*F(1)-E(1)*F(2)*A18 

A83=A82*A82 

A84=A81+A83 

A22=DSQRT(A84) 

A23=(A20*A17)/A22 

A24=l.  0D0-A23*A23 

A90=A22/A21 

A25=DATAN(A90) 

IF(A25)  80,85,85 
*"  80  A25=(PI/2.0D0)+A25 
85  A26=A25*A25 

A27=l.  0D0/A22 

A91=A21/A22 

A29=A24*A24 

A30=(A50*A25)+A19*(A51*A22-A52*A26*A27) 

A31=A24*(A53*A25+A53*A22*A21+A52*A26*A91) 

A32=A19*A19*(-A52*A21*A22) 

A33=A29*(A54*A25+A54*A22*A21-A52*A26*A91-A55*A22*A21**3) 

A34=A19*A24*(A52*A26*A27+A52*A22*A21*A21) 

S=(A30+A31+A32+A33+A34)*B 
C 

C     COMPUTATION  OF  DIFFERENCE  OF  LONGITUDE  ON  THE  REDUCED  SPHERE 
C 

A36=(A51*A25)+A19*(-A52*A22-.A56*A26*A27) 

A37=A24*(-A57*A25+A58*A22*A21+A56*A26*A91) 

A38=(A36+A37)*A23+A16 


c 

c 
c 

COMPUTATION  OF  GEODETIC 

FORWARD  AZIMUTH 

A39=DSIN(A38) 

A40=DC0S(A38) 

A41=(E(2)*F(1)-A40*E(1)- 

'•-F(2))/(A39*F(2) 

IF(A41)  95,90,95 

90 

A41=0. 5D0-7 

95 

A42=l. 0D0/A41 
A43=DATAN(A42) 
IF(C35)  120,100,100 

100 

IF(C35-PI)  105,115,115 

105 

IF(A41)  110,140,140 

110 

A43=PI+A43 
GO  TO  140 

115 

IF(A41)  130,135,135 

120 

IF(C35+PI)  105,105,125 

125 

IF(A41)  130,135,135 

130 

A43=PI-A43 
GO  TO  140 

135 

A43=2.  0D0*PI-A43 

140 

A43=A43+PI 
A43=A43-2. 0D0*PI 
IF(A43)  145,150,150 

145 

A43=A43+2.  0D0*PI 

150 

FAZ=A43/ARC1 
FAZ=FAZ+648000.  0D0 

IF(FAZ-1296000.0D0)  160 

,155,155 
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155  FAZ=FAZ-1296000.0D0 
C 

C     COMPUTATION  OF  BACK  AZIMUTH  NOT  INCLUDED  IN  THIS  VERSION 
C 

160  BAZ=0. ODO 
RETURN 
C 

C     PRINT  ROUTINE  FOR  INVALID  DATA 
C 

255  PRINT  270 

CALL  ANGLE(3,ALAT1,ISGN1,LAT1D,LAT1M,ALAT1S) 
CALL  ANGLE(3,AL0N1,ISGN2,L0N1D,L0N1M,AL0N1S) 
CALL  ANGLE ( 3 , ALAT2 , I SGN3 , LAT2D , LAT2M , ALAT2 S ) 
CALL  ANGLE ( 3 , ALON2 , I SGN4 , LON2D , LON2M , ALON2S ) 
265  PRINT  280,ISGN1,LAT1D,LAT1M,ALAT1S,ISGN2,LON1D,LON1M,ALON1S,ISGN3, 
*LAT2D , LAT2M , ALAT2S , I SGN4 , L0N2D , LON2M , ALON2S 
RETURN 
C 
C     FORMAT  STATEMENTS 

270  F0RMAT(1H1,51X,' SUBROUTINE  SODIN  INVALID  DATA*/) 
280  FORMAT( IX, 4(2X,A1, 213, F7. 3), 2X, 12, 2X,' IDENTICAL  STATIONS') 
END 


SUBROUTINE  ANGLE  (N,ARGS , ISIGN, IDEG,MIN,SEC) 


IMPLICIT  REAL*8(A-H,0-Z) 

DIMENSION  ASEC(6) 

DATA  ASEC/59. 45D0,59. 945D0,59. 9945D0,59. 99945D0,59. 999945D0, 
*59. 9999945D0/ 

DATA  MINUS/ 1H-/ , IBLANK/1H  / 
C 

C     ANGLE  CONVERSION 
C 

IF  (ARGS)  10,5,5 
5    ISIGN=IBLANK 

BRGS=ARGS 

GO  TO  15 
10   ISIGN=MINUS 

BRGS=-ARGS 
15   IDEG=BRGS/3600.0D0 

ARGT=IDEG*3600 

AMIN=BRGS-ARGT 

MIN=AMIN/60. ODO 

AMINT=MIN*60 

SEC=AMIN-AMINT 

M=N+1 

IF  (SEC-ASEC(M))  30,20,20 
20   SEC=0. ODO 

MIN=MIN+1 

IF  (MIN-60)  30,25,25 
25   MIN=0 

IDEG=IDEG+1 
30   RETURN 

END 
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APPENDIX  D.     TABLE  FOR  STANDARD  DEVIATIONS  AND 
CO  VARIANCES  OF  VELOCITY  COMPONENTS  FOR  MINI  RANGER 


TIME  IN- 
TERVAL 

ov 

<>u 

Vy.u 

3  MIN 

1.3  cm/sec 

1.3  cm/sec 

0.1234  cm2/  sec2 

6  MIN 

0.64  cmsec 

0.64  cm/sec 

0.0386  cm2/  sec2 

12  MIN 

0.32  cm,  sec 

0.32  cm/sec 

0.0116  cm2/ sec2 

20  MIN 

0.19  cm  sec 

0.19  cm/ sec 

10-1  x  0.0486c//?2/  sec2 

40  MIN 

10"3x  0.98c/??/ sec 

10~3x  0.97 cm/  sec 

10-2x  0.17c/??2/ sec2 

60  MIN 

10~3x  0.67c/??/  sec 

10-3  x  0.66cm/  sec 

10"2x  0.1080c/??2/  sec2 
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APPENDIX  E.  AVERAGE  CURRENT  PROFILES 
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Figure  17.      AVERAGE  CURRENT  PROFILE  FROM  POINT  1  TO  POINT  3 
USING  NAVIGATION  DATA  FROM  MINI  RANGER 
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Figure   18.       AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO  POINT  2 
USING  NAVIGATION  DATA  FROM  MINI  RANGER 
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Figure  19.       AVERAGE  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  3 
USING  NAVIGATION  DATA  FROM  MINI  RANGER 
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Figure  20.       AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO  POINT  4 
USING  NAVIGATION  DATA  FROM  MINI  RANGER 
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Figure  21.      AVERAGE  CURRENT  PROFILE  FROM  POINT  1  TO  POINT  3 
USING  NAVIGATION  DATA  FROM  LORAN  (TD) 
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Figure  22.      AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO  POINT  2 
USING  NAVIGATION  DATA  FROM  LORAN  (TD) 
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Figure  23.      AVERAGE  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  3 
USING  NAVIGATION  DATA  FROM  LORAN  (TD) 
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Figure  24.       AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO  POINT  4 
USING  NAVIGATION  DATA  FROM  LORAN  (TD) 
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Figure  25.      AVERAGE  CURRENT  PROFILE  FROM  POINT  1  TO  POINT  3 
USING  NAVIGATION  DATA  FROM  LORAN  (DISPLAY) 
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Figure  26.      AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO  POINT  2 
USING  NAVIGATION  DATA  FROM  LORAN  (DISPLAY) 
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Figure  27.       AVERAGE  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  3 
USING  NAVIGATION  DATA  FROM  LORAN  (DISPLAY) 
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Figure  28.      AVERAGE  CURRENT  PROFILE  FROM  POINT  3  TO  POINT  4 
USING  NAVIGATION  DATA  FROM  LORAN  (DISPLAY) 
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APPENDIX  F.     CURRENT  PROFILES  FOR  SPECIFIC  TIMES 


Figure  29.  CURRENT  PROFILE  FROM  POINT  1  TO  POINT  2  USING 
NAVIGATION  DATA  FROM  MINI  RANGER  FOR  A  SPECIFIC 
TIME  (14  33  11) 
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Figure  30.  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  1  USING 
NAVIGATION  DATA  FROM  MINI  RANGER  FOR  A  SPECIFIC 
TIME  (15  06  10) 
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Figure  31.  CURRENT  PROFILE  FROM  POINT  1  TO  POINT  2  USING 
NAVIGATION  DATA  FROM  LORAN  (TD)  FOR  A  SPECIFIC 
TIME  (14  33  11) 
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Figure  32.  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  1  USING 
NAVIGATION  DATA  FROM  LORAN  (TD)  FOR  A  SPECIFIC 
TIME  (15  06  10) 
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Figure  33.  CURRENT  PROFILE  FROM  POINT  1  TO  POINT  2  USING 
NAVIGATION  DATA  FROM  LORAN  (DISPLAY7)  FOR  A  SPE- 
CIFIC TIME  (14  33  11) 
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Figure  34.  CURRENT  PROFILE  FROM  POINT  2  TO  POINT  1  USING 
NAVIGATION  DATA  FROM  LORAN  (DISPLAY)  FOR  A  SPE- 
CIFIC TIME  (15  06  10) 
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